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Abstract 

The goal of the Fully Online Datacenter of Ultraviolet Emissions (FONDUE) Work- 
ing Team of the International Space Science Institute (ISSI) in Bern, Switzerland, was to 
establish a common calibration of various UV and EUV heliospheric observations, both 
spectroscopic and photometric. Realization of this goal required a credible and up-to-date 
model of spatial distribution of neutral interstellar hydrogen in the heliosphere, and to 
that end, a credible model of the radiation pressure and ionization processes was needed. 
This chapter describes the latter part of the project: the solar factors responsible for shap- 
ing the distribution of neutral interstellar H in the heliosphere. Presented are the solar 
Lyman-alpha flux and the question of solar Lyman-alpha resonant radiation pressure force 
acting on neutral H atoms in the heliosphere, solar EUV radiation and the process of pho- 
toionization of heliospheric hydrogen, and their evolution in time and the still hypothetical 
variation with heliolatitude. Further, solar wind and its evolution with solar activity is 
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presented, mostly in the context of the charge exchange ionization of heliospheric neutral 
hydrogen, and in the context of dynamic pressure variations. Also the question of electron 
ionization and its variation with time, heliolatitude, and solar distance is presented. After 
a review of the state of the art in all of those topics, we proceed to present an interim 
model of solar wind and the other solar factors based on up-to-date in situ and remote 
sensing observations of solar wind. This model was used by Izmodenov et al. to calculate 
the distribution of heliospheric hydrogen, which in turn was basis for intercalibration of 
the heliopsheric UV and EUV measurements discussed in the other chapters of the book. 
Results of this joint effort will further be utilised to improve on the model of solar wind 
evolution, which will be an invaluable asset in all heliospheric measurements, including, 
among others, the observations of Energetic Neutral Atoms by the Interstellar Boundary 
Explorer (IBEX). 



2.1 Brief description of the physics of neutral interstellar 
gas in the inner heliosphere 

Distribution of neutral interstellar hydrogen and ultraviolet radiation in the inner he- 
liosphere are closely interrelated. Absolute calibration of observations of the heliospheric 
backscatter Lyman-alpha glow requires both the knowledge of the well-calibrated solar 
EUV output and of other solar forcing factors, mainly the solar wind. The role of those 
factors and their variabilities in shaping the distribution of neutral interstellar hydrogen 
can be apprehended from modeling papers cited in the remaining portion of this section. 

If one takes an inflow with a finite velocity of a fully neutral gas at a temperature 
Tea far away from the Sun, then when the ionization rate /3 at r and the effective force 
F (r) acting on the atoms are spherically symmetric and time-independent, the distribution 
function of the gas at a distance r from the Sun will feature axial symmetry around the 
inflow direction and be given by the equation: 

v • V r f (v, r) + ^ • V v / (v, r) = -pf (v, r) (2.1) 
m H 

where V v is the gradient operation in the x-space, r and v are, respectively, position and 
velocity vectors of the gas cell element, and m# is hydrogen atom mass. Together with 
the assumption that the gas "at infinity" (in practice: a few hundreds of AU from the Sun) 
is homogeneou s and Maxwellian and t hat the gas is collisionless at the distances relevant 
for the model uzmodenov et al. , 2000] this is the basis for the classical hot model of the 



distribution of neutra l interstellar gas in the he liosphere I r/ioma5lll978t[Fh/»lll978lll979l: 



Wu and Judee . 19791: Lallement et al l | l985all . Based on Liouville theorem, solution of 



this equation | Danby and Camni 1957ll forthe distribution function / (v, r, t) for time t 



location r and velocity v can be expressed as: 

/ (v, r, t) = /«, (Vco (v, r) , r ra (v, r)) W (v, r, t) , (2.2) 

where W (v, r, f) is the survival probability of an atom that at the time t and location r 
has velocity v and which at a distant location had velocity v^. For now, t is only a 
formal parameter here. The probability of existence of such an atom in the distant region 
of the heliosphere is given by the distribution function (v M , r ra ), and the link between 
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the local velocity and position vectors v (/), r (/) and the corresponding vectors in the so- 
called source region of the atoms can be obtained from solution of the equation of motion 
of hydrogen atoms in the heliosphere: 

F(r,t,v r) = - GmM(l -/ (V " t))r - (2.3) 

where F is the total force acting on the atom with a mass m, G is gravitational constant, 
M solar mass, v r radial velocity of the atom at time t, and ju is the ratio of solar resonant 
radiation pressure force to the solar gravity, which will be discussed further on in this 
chapter. 

As can be seen from this description, the distribution of neutral interstellar hydrogen 
in the inner heliosphere is determined on one hand by the dynamical influence of the Sun 
by a balance between its gravity and radiation pressure forces and on the other hand by 
the ionization losses. Both will be extensively discussed in the following parts of this 
chapter. Here we only say that the ionization processes include charge exchange between 
the incoming neutral atoms and solar wind protons, ionization by impact of solar wind 
electrons, and ionization by the solar EUV radiation. 

The hot model, however handy it is being almost analytical (in fact, numerical calcula- 
tions are needed only when integrating the local distribution function to yield its moments, 
like density and mean velocity), is not perfect because none of its assumptions is ful- 
filled. First, the interstellar gas in the Local Cloud is not fully neutral and the interaction 
of its ionized component with the plasma of the solar wind creates a boundary region of 
the heliosphere (the heliospheric interface), which begins at the termination shock of the 
solar wind, where the solar wind becomes subsonic and eventually turns back at the he- 
liopause, which can be approximated as a thin layer separating the solar wind plasma from 
the interstellar plasma. Beyond the heliopause there is the outer heliosheath, where the 
pristine neutral interstellar gas becomes altered due to charge exchange interactions with 
the protons from the piled and heated up interstellar plasma. History of the de velopment 
of modeling of this region of the heliosphere can be found in Bamnov [2006b] and mod- 
ern views on this topic have b een recently reviewed by Fahr [2004]; B aranow 1 2006a]; 



Izmodenov and Bara nov [ 2006] and by Izmodenov et aL in this book. 

Quantitative interpretation of the helioglow measurements required improvements in 
the cla ssical hot model, which were realized quite early in the history of heliospheric re- 



search. Lal lement et al. [1985b] allowed for latitudinal modulation of the charge exchange 



rate, approximating it with a one-parameter formula 1 - A sin 2 <p, which enabled them to 
vary the equator-to-pole ratio of the ionization rates, but required keeping the width and 

range of the equatorial region of enhanced ionization fixed; 

A different extension of the hot model was proposed bv lRucinski and Fahr 1 198911 199 ill 



who pointed out that the rate of electron-impact ionization does not change with r" 1 , even 
though its effects on the distribution of neutral interstellar hydrogen in the heliosphere 
are noticeable only within a few AU from the Sun, where its density is already strongly 
reduc ed. Therefore t his aspect of the heli ospheric physics has been left out afterwards, 
until Bzowski | 2008 1 ; \Bzowski et all [ 2008ll reintroduced it in a refined, latitude-dependent 



manner (see also further down in this chapter). 

The next generation of heliospheric models abandoned the assumption of invariability 
of the solar radiation pressure forc e and of the ionizati on rate. The first, though very much 



simplified model was proposed bv \Kvrdla et al.\ 11994], followed bv \Ruciriski and Bz.owski 
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1 1995a |; Bzowski and Rucinski 1 1995al IbTl: \Bzowski et al. 1 1997 1. who studied variations 
in density, bulk velocity, and temperature of neutral interstellar hydrogen near the Sun, as 
well as variations in the helioglow intensity assuming (lacking sufficient data at that time) 
an analytic model of evolution of the radiation pressure and ionization rate during the solar 
cycle. 



Another modification wa s introduced b y Scherer et al. [ jl999ll , w ho following Osterba rt and Fahn 
1 19921: Ba ranov et al. ; Baranov and Malama 1993] [see also Malama et al. 2006] 

realized that the charge exchange processes at the boundary layer of the heliosphere cre- 
ate a new, so-called secondary collisionless population of neutral H atoms, which can 
be taken into account by approximating the distribution function in Eq. l2.2l bv a sum 
of two two-dimensional Maxwellian functions with anisotropic temperatures. The val- 
ues of temperatures, densities and bulk velocities in these Maxwellian functions, how- 
ever, had to be adopted as some functions of the angular separation of the point in 
Eq. |2^2| fro m the upwind direc t ion. S uch an approach was later expanded and improved 
by \Katushkina and IzmodenovX 1201011 . In this approach, however, the parameters of the 
distribution function /«, must be obtained from a n external model, such a s the Moscow 
Monte Carlo model of the heliosphere [see, e.g., \lzmodenov et al. , 2009]. Along with 



the two-populations non-Gaussian model, an approximation of the radiation pressure and 
ionization rate by spherically symmetric series of sines and cosines fitted to actual mea- 
surements of the ionization rate a nd radiation p r essure in t he ecliptic was a dded. These 
approximations were described by Sch erer et al. |1999] snd \Bzowsk i |2001a, b]. 



The most recent developmen t in the modeling was addition by uarnopols ki [2007]; 
Tarnopolski and Bzowski 1 2008a] of the radiation pressure force being function of radial 
velocity of an atom relative to the Sun; this effect will be discussed further in this chapter. 



The list of modifications to the classical hot model presented above is also a (perhaps 
not full) list of effects that need to be taken into account at the solar side to facilitate inter- 
calibration of the measurements of the helioglow with other UV observations in space. 
Apart from the heliospheric side, there is also physics of the heliospheric inter face and the 
conditions in the Local Interstellar Cloud ([see, e.g., \Frisch et al. . I2009ll201ll for review] 
that must be taken into account, which, however, are beyond the scope of this chapter. 



From the above description it is clear that accurate modeling of neutral interstellar hy- 
drogen in the inner heliosphere requires accurate knowledge of the factors contributing to 
the ionization and radiation pressure. For the radiation pressure force one needs to know 
the evolution of the shape of the solar Lyman-alpha spectral line in time. The ioniza- 
tion processes include charge exchange with solar wind protons (and perhaps also alpha 
particles), ionization by electron impact, and ionization by solar EUV radiation. Recom- 
bination is not important here because of two reasons: (1) its rate is small as compared 
with the ionization rate [ Wachowicz, 2006], and (2) the recombination products maintain 
their velocity, which is equal to the solar wind velocity, so they do not contribute to the 
populations of heliospheric atoms that are capable of scattering the solar EUV radiation 
responsible for the helioglow. 



2.2. Radiation pressure 
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Figure 2.1: Wavelength- and disk-integrated solar Lyman-alpha flux from the Laboratory 
for Atmospheric and Space Physics (LASP), referred to as the total Lyman-alpha flux 7 tot . 
The daily time series is a composite of actual measurements from various experiments re- 
scaled to the common calibration of UARS/SOLSTICE, with the gaps filled by the proxies. 
Color codes: gray: F10.7, orange: AE-E, yellow: SME, green: MgII c / w proxy, aquamarine: 
UARS/SO LSTICE version 1 8, light blue: TIMED/SEE, dark blue: SORCE/SOLSTICE. 
Based on Woods et al. 



2.2 Radiation pressure 

2.2.1 Evolution of the total solar flux in Lyman-alpha in the ecliptic 

The radiation pressure force acting on neutral interstellar H atoms in the heliosphere 
is proportional to the total flux in the solar Lyman-alpha spectral line, which is defined as 
the spectral flux integrated ove r a 1 nm interval fro m 120 to 121 nm. It has been measured 
since the middle of 1970s \see \W00ds et al. , 2000, for the history of measurements]. De- 
spite all efforts, while precision of the measurements was good, the problem of absolute 
calibration, prone to changes with time, has affected the accuracy from the very beginning 
until present and a measure of progress in this field is the reduction of discrepancies from a 
factor of 4 in the 1970s to ~ 15% nowadays. The total composite Lyman-alpha time series, 
available from the Laboratory for Atmospheric and Space Physics (LASP) at the Univer- 
sity of Colorado in Boulder, Co, is scale d to the absolute calibration of UARS/SOLSTICE 



I Woods et al. , 19961: rWfootfc et fl/.ll2000ll . The cadence of measurements is presently 1 day 



and the inevitable gaps are usually filled using a hierarchy of proxies, as illustrated in 
Fig.O 

The most widely used is th e proxy based on the solar radio flux at the 10.7 cm wave 



length, the so-called F 10.7 flux \ Covinston , 19471: TbppiniA. 1 1987ll. Another frequently used 
is the Mg II core-to- wing (MgII c / w ) index IHeath and SchlesineeA\l98(i . \Viereck andPueal 



1999h . 



The use of proxies raises the question of credibility of the results \Flo vd et q/.L l2002. 2005] 
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Figure 2.2: Rates of change of the total Lyman-alpha flux 7 tot . Shown are the rates of 
change per unit time ^f- = ^afeiz^aw G f t ne daily (upper panel) and monthly (lower 
panel) composite Lyman-alpha flux shown in Fig. 12.11 with identical color coding. 



The solar EUV radiation is variable on many time scales, from hours to more than solar 
cycle length. Proxies generally follow the variability of the quantity being approximated, 
but not precisely and not at all time scales. In particular, even though the correlation co- 
efficients for daily values may exceed 0.9, at the shortest time scales the agreement can 
be problematic because - depending on the location of the source region of a given proxy 
quantity - they can show the effect of limb brightening or l imb darkening different from 
the quantity they are used to approximate IFlovd et all 1200511 : an illustration can be found 
in the upper panel of Fig. 12.21 In consequence, the quality of the approximation at these 
short scales is reduced. 

On the other hand, the rates averaged by Carrington rotation (alternatively referred to 
as "monthly" in this text) seem to be reproduced quite well, as judged from the behavior 
of the monthly change rat es shown in the lower panel of this figure and as suggested by 



ig' 
12 



Dud ok de Wit et al. [2009]. For modeling of the distribution of neutral interstellar gas in 
the heliosphere the Carrington period is anyway the finest time scale presently in use, so 
the proxy imperfections are not a big problem in this respect. 



2.2.2 Variation of the Lyman-alpha flux with heliolatitude 



T he solar disk-averaged Lyman-alpha flux is made of at least three components mmblard et al. 
2008]: from the quiet Sun, from the coolest part of the chromosphere, and from the 
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hot lower coro na. The inhomoge neous distribution in heliolatitude of active regions was 
pointed out by Cook et al. 1 1980tl . who construct ed a two-c o mpon ent latitude-dependent 
model of disk-averaged solar UV irradiance. \Cook et al.\ M 1 93 in considered the solar 



Lyman-alpha emission and suggested that the ratio of the disk-integrated solar flux at the 
pole to the flux at the equator should be about a^ xa = 0.8 during solar minimum, when 
the active regions are distributed in latitudinal bands. Thes e suggest i ons w ere supported 
by direct solar minimum observations of the solar c orona b y Auchere J2005h . Such a ratio 
was suggested as valid also for solar maximum b\ \Prvor et all [1992], based on indirect 
evidence from observations of the Lyman-alpha helioglow. Thus it seems that the latitu- 
dinal anisotropy of the Lyman-alpha flux does not change substanti ally during the solar 



cycle, although this conclusion certainly needs further verification. Wzow ski [2008] sug- 
gested that the latitude dependence of the disk-integrated solar Lyman-alpha flux may be 
approximated by the formula 



: (<P) = hot (0) ^fl L ya SHI 2 (<p) + COS 2 (0) 



(2.4) 



where <p is heliolatitude, flLya the "flattening" factor and / tot (0) is a equatorial Lyman-alpha 
flux. 

Heliolatitude variation of the disk-integrated Lyman-alpha line profile is, to our knowl- 
edge, un explored. On one hand, one should expect some form of variation because, as 
shown by Tian et al. l ll2009aH bLlcTl. the line profile depends on the features on the solar disk 
that are being observed and the latitude distribution of these features is inhomogeneous 
and varies during solar cycle. On the other hand, the apparent lack of strong variability of 
the disk-integrated flux might suggest that the spectral variation with helio latitude is mild. 



Recen t investigation of the variation of solar spectrum with heliolatitude bv \Kiselman et al. 



1 201 1] seems to bring a negative result (i.e., no variation). Should the disk-integrated spec- 
tral flux indeed vary with heliolatitude, this would potentially have consequences both for 
the photoionization rate and the radiation pressure force. 



2.2.3 Mechanism of radiation pressure 

Mechanism of resonant interaction of a H atom with solar radiation, which leads to 
the force of resonant radi ation pressure acting on ne utral H atoms in the heliosphere, was 



extensively discussed by Bras ken and Kvrola\ II 1 99811 . In brief, the probability f aD a (A) that 



a hydrogen atom in the ground state, whose base wavelength is Ao, absorbs an incoming 
photon at wavelength A is equal to: 

/ abs (A) = r * r2 (2.5) 



+ 



where Yr is the energetic width of the second orbital of the atom, corresponding to about 
±25 m s around the base wavelength of the Lyman-alpha transition. The atom whose 
radial velocity relative to the Sun is will absorb photons from the very center of the solar 
line, but if its radial velocity v r is non-zero, then due to Doppler effect it will be tuned to 
a different portion of the solar line profile, namely to the wavelength A - Aq (1 - v r fc). 
Within about 10~ 15 s after absorption, the atom will re-emit the Lyman-alpha photon at an 
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angle u to the impact direction with the scatter-angle probability p {to) described by: 



p(oj) 



cos (lo) 1 1 
4 + L2 



(2.6) 



Hence a resonant interaction of the atom with a suitable photon results in a change of 
atomic momentum at the moment of absorption by Ap = hv = c h/A in the antisolar di- 
rection, followed after a time of 1CT 15 s by another momentum change in the direction 
described by Eq. 12.61 However, typical interactions frequencies at 1 AU, which are pro- 
portional to the solar spectral flux, are on the order of 1/500 Hz I OuemeraiX l2006ll . so 
statistically at time scales which are shorter than the time scales of a change in atomic 
velocity relative to the Sun the only net effect of the interaction of the atom with solar 
radiation is the antisolar momentum change. Since the in t erpla netary medium is almost 
optically thin within a few AU from the Sun Wuemerais . l200dl . the solar spectral flux 
scales with inverse square of solar distance and consequently also the solar radiation pres- 
sure force scales with solar distance as 1/r 2 , which gives effect of compensation of solar 
gravity by a factor p. Hence the effective solar force acting on the atom is conveniently 
expressed as the factor p (I tot (t) , v r ) of solar gravity. It is proportional to the spectral flux 
Fa corresponding to the Doppler-shifted wavelength A = Aq (1 + V r /c), resulting from the 
instantaneous radial speed v r of the atom relative to the Sun. Since the spectral flux varies 
with time, effectively the p factor is a function of radial velocity and time, as expressed in 
Ea. 1231 



2.2.4 Profile of solar Lyman-alpha line and the solar radiation pres- 
sure factor jx 

Measu rements of the solar Lyman-alpha line profile, although scarce, date back to 



the 1970s | Vidal-Madiar 1975; A rtzner et all [1978; Bon net et ail 11978c \Lemaire et al. 



19781: \Woods et al.[ Il995ll . They were performed from within the Earth exosphere and 



hence suffered from the absorption by geocoronal hydrogen in the spectral region most 
relevant for the helioglow. Only after launch of SOHO, which orbits at the LI Lagrange 
point, an u nobscured view on the full spectral range of the d i sk- integrated solar lin e was 
obtained bv marren et al\ M1998al ; \Lemaire et al.\ 111998 , 2002 ] ; Lemaire et al.\ M2005I1 , who 
used a novel, ingenuous observations technique proposed by J.-L. Bertaux. Owing to the 
latter measurements it was possible to study the evolution of the profile at all solar activity 
phases from minimum to maximum. 

The solar Lyman-alpha li ne features a s elf-reversed shape that used to b e approximated 



1997; Scherer et al. 



by tw o Gaussian functions WahA 1 1979c \C habrillat and Kockarts, _ 
2000]. Recently, Tarnovolski 1 2007 |] : Uarnopolski and Bzowski 1 2008b] showed that the 
measurements bv \Lemaire et al.\ [2002] can all be fit by a three-Gaussian model parame- 
terized by the disk-integrated flux: 

H (v„ /tot (0) = A [1 + B/ tot (?)] exp (-Cv 2 r ) [l + Dexp [Fv r - Gv 2 r ) + i/exp {-Pv r - Qvfj] 

(2.7) 

with the following parameters: 



A = 2.4543 x 10~ 9 , 
D = 0.73879, 
H = 0.47817, 



B = 4.5694 x 10~ 4 , 
F = 4.0396 x 10~ 2 , 
P = 4.6841 x 10~ 2 , 



C = 3.8312 x 10~ 5 
G = 3.5135 x 10~ 4 
Q = 3.3373 x 10~ 4 
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Compensation factor fi as function of radial velocity 
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Figure 2.3: Factor fi of solar radiation pressure force to solar gravity based on the model 
specified in Eq. 12.71 shown as a function of radial velocity of a H atom relative to the Sun 
for the total solar flux values corresponding to the minimum (red) and maximum (blue) of 
solar activity. Thick lines indicate the spectral region ±30 km s around Doppler shift, 
relevant for the neutral interstellar hydrogen gas in the heliosphere. 



The accuracy of the fit is similar to the accuracy of the measurements, assessed to be 
~ 10%. 

With this formula, one can calculate the fi factor for an arbitrary radial velocity v r pro- 
viding that the total solar Lyman-alpha flux 7 tot is known. The dependence of the fi factor 
on radial velocity for the total flux values repre sentative for solar minimum and maximum 
conditions are shown in Fig. 12.31 adapted from lTarnopolski and Bzowski 1 2008bll . 

The spectral region of the solar Lyman-alpha line most relevant for modeling of the 
heliospheric Lyman-alpha glow in the wavelength band straddling the central wavelength 
by approximately ±30 km s" 1 . The spectral flux at the line center is closely correlated 
with the line-integrated flux h n \Vidal-Madiar and PhissamaA I "l980ll . This is beneficial to 
the modeling of the helioglow because it permits to easily calculate the fi factor based on 
measurements of the solar line-integrated Lyman-alpha flux. 



2.2.5 Approximate values of the /u factor 

Since modeling of neutral interstellar hydrogen with the use of a model that takes the 
full solar line profile is computationally demanding, the approach where the fi factor is 
treated as independent of radial velocity are still widely used. To that purpose a formula 
to translate the line-integrated flux into the fi factor is needed. In the past, this issue was 
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Solar gravity compensation factor fi 
as function of total flux 
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Figure 2.4: Solar radiation pres sure factor n shown a s a function of the total flux in the 
Lyman-alpha line calculated by \Emerich et al. 1 2005 1 from Eq. 12.91 which connects the 
total flux with the spectral flux precisely at line center (blue), and b\ \Bzowski et al\ [2008] 
from Eq. 12.101 (red), where correlation between the total solar flux and spectral flux is 
adopted from averaged spectral flux over +30 km s _I around the line center. 



addressed just by taking 



//(/ tot ) = 3.0303 x 10 10 a/t 



(2.8) 



where 7 to t is the disk- and line-integrated solar Lyman-alpha flux and a is a constant usu- 
ally adopted between 0.85 and 1. Wit h increasing accuracy of measurements, more so- 
phisticated formulae became available. \Emerich et al.\ 11200511 fitted the following relation 
between the spectral flux at line centre to the total flux: 



10 12 cm 2 s 'nm 



= °- 64 ( inii / T 2 r' ±o - 08 

\10 n s 'cm 2 / 



(2.9) 



while Wzowski et all 11200811 found a linear relation between the spectral flux averaged by 
the spectral range ±30 km s about the line center and the line- and disk-integrated flux: 



fJ- (hot) 



3.4731Q- 12 /,, 



0.287 



(2.10) 



A comparison of the two predicted fi values as a function of the solar total Lyman-alpha 
flux is shown in Fig. 12. 41 and suggests that if one decides not to use the model of radiation 
pressure force dependent on the radial velocity of the atom, then calculation of a good 
effective fi factor is not a straightforward task. Use of the model with n being a function of 
v r is ob ligatory when one wants to model interstellar deuterium I Tarnopolski and Bz.owski , 
20083. 



The // values returned by the formulae bv \Emerich et al.\ B2005I1 and uarnopolski and Bz.owski 
1 2008b] calculated from the monthly values of the LASP composite Lyman-alpha flux 
shown in Fig. 12. II are presented in Fig. 12.51 Differences between the /j values obtained 
from equations d2.9l l and ( 12.101 1 are on the order of 10%, i.e. on the order of the uncer- 
tainty of the total flux. 



2.3. Ionization processes 
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Figure 2.5: Comparison of the solar radiation pressure factors [i that approximate the 
compensation of solar gravity by the resonant radiation pressure force acting on neutral H 
atoms. Gray dots are daily values from the composite Lyman-alpha flux series presented 
in Fig. 12. 11 blue line is the monthly (Carrington rotation averaged) val ues of the flux scaled 
using Eq. 12.81 green line are monthly values obtained from Eq. 12.91 \Emerich et aZ [ | 2005ll 
and red line represents monthly values obtained from Eq. 12. 101 1 Bz.o wski et a/.[l2008ll . 

2.3 Ionization processes 

The three main ionization processes of neutral interstellar hydrogen atoms in the he- 
liosphere are the following: 

charge exchange with solar wind charged particles (mostly protons): 

H + p — > ppui + Hena 
photoionization by solar EUV radiation: 

H + v — > ppui + e 
and ionization by impact of solar wind electrons: 

H + e — > ppui + 2e 

The reaction product is always a pickup proton ppui and either an energetic neutral atom 
Hena or one or two electrons. There are no recombination processes that would inject 
neutralized protons to the region of phase space occupied by the population of interstellar 
neutral atoms. 

As a result of the charge exchange reaction, a solar wind proton captures the electron 
from the neutral interstellar H atom and becomes an energetic neutral atom (ENA), but 
maintains its momentum and thus does not enter the interstellar population. In this respect, 
even though one neutral H atom is replaced by another, such a reaction is still a loss process 
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despite the total number of H atoms in the system does not change - but it does change in 
the population under study. 

The newly-created protons regardless of the reaction they orig inate from, are picked up 



1 pr 

by the solar wind flow IFq/tH,! 19731:1 Vasvliunas andSiscoeW 19761], cre ating a distinct popu 



lation that can be measured \Mbbius et al. , 1985; Gloeckler et all Il993h and analyzed [e.g. 



Gloeckler and G eiss. 2001; Gl oeckler et ali 12004: Bzo wski et all 1200811 . but this problem 
is outside the scope of this text. The solar wind is not depleted in protons due to the charge 
exchange reaction since eventually a solar wind core proton is replaced with a pickup pro- 
ton - but the concentration of protons per unit volume is not changed even though the 
distribution function of solar wind protons is modified. 

The two remaining ionization reactions cause genuine losses, the H atom enters the 
reaction and is not replaced with another at a different velocity: eventually, a proton- 
electron pair is created. 



2.3.1 Charge exchange 
General formula 

The process of resonant charge exchange between H atoms and protons is of cru- 
cial importance for the physics of the heliosphere. It contributes to the pressure balance 
between the solar wind and interstellar gas and makes possible energy and momentum 
transport between solar wind and interstellar gas across the heliopause. Charge exchange 
losses of neutral interstellar gas in the supersonic solar wind, which will be discussed in 
this section, are only a small piece of a large picture. 

The rate of charge exchange between neutral H atoms and solar wind protons can be 
regarded as probability of a charge exchange act per unit time in a given location in space. 
For a H atom traveling with velocity v# and a local proton distribution function f p (v p ), 
where v p is the velocity vector of an individual proton, the rate of charge exchange can be 
calculated from the formula: 

Pcx = J o-cx(\vH-v p \)\v H -Vp\fp(yp)dv p , (2.11) 

where ctqx (\vh - V P U is the reaction cross section and v H -v p = v re i is the relative velocity 
between the H atom and an individual proton. The integration covers the entire proton 
velocity space. This formula can be put into an equivalent form: 

Pcx = J ocx (|v re i|) |v re i| f p (v H - v te i) dv re i (2.12) 

Depending on the underlying plasma regime and on the velocity of the H atom, various 
simplifications can be made. When the kinetic spread of the plasma ui 4 , is small compared 
with the plasma flow velocity Vsw: ut, p = ^2kT p /m p <K vsw = fv P f (vpj dv p /n p , we can 
approximate the proton distribution function by delta-function centered at the solar wind 
speed and the formula for charge exchange rate simplifies to: 

PCX = OCX (Vrel) n p v leh (2. 13) 

where n p is the local proton density and v re i becomes v re i = |v# - v$wl This is the case 
for the scenario of Hena in the supersonic solar wind for the atoms that travel at vh ~ 
50 km s _1 or faster. 
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For vh <S vsw> as in the case of the thermal interstellar H populations in the supersonic 
solar wind, we have v re i = v p . Then the rate of charge exchange between H atoms and solar 
wind protons will be given by 

jScx = 0-cx (vsw) « P Vsw (2.14) 

This is not the case, however, for Hena which travel at vh >50 km s _1 , which is not 
negligible compared with the typical solar wind speed of vsw ~440 km s _1 , in the subsonic 
solar wind (inner heliosheath), and within the Baranov wall (outer heliosheath) thermal 
speeds of local plasmas are comparable to specific velocities of H atoms and the full form 
of the charge exchange formula must be applied. 



Charge exchange cross section 



The collision speed range most relevant for the heliospheric physics is from ~ 1 km s 
to ~ 1000 km s _1 , which is equivalent to 0.005 eV and 5.2 keV. Relative speed be- 
tween interstellar neutral H atoms and protons in the supersonic solar wind range from 
~ 300 km s- 1 to ~ 1000 kms 4 . 

A detailed discussio n of the cha r ge exc hange process and of the cross sections for this 
reaction can be found m Wahr et al\ M2007I1 and will not be repeated here. For the purpose 
of this work it is important to point out that there were basi c ally 4 cross section formu - 



lae used in the heliospheric physics: fro m Fite et al [ 1962]; Mahe r and Tinslew B1977I1 : 



Barnett et al. 



Fite et al. 



199Q]; \Lindsav and StebbinsT l 2005 



5 7 ', 



jl962ll and Maher and TinsleM | 1977H proposed to approximate the charge 
exchange cross section as a function of relative velocity v re i between the colliding partners 
by a formula 

ocx (v re i) = (a + b In v re i) 2 (2. 15) 

The v alidity of the formu la by Fite et al. i cm 

was claimed to be between 20 and 2000 eV 
and of Maher and Tinsley 1 1977] between 0.005 and 1 keV. 



Barnett et al.\ 11990] prop osed and involv ed fit in a form of Chebyshev polynomials 



valid i n a broad energy range. \Bzowski [2001b] approximated the data used bv uiarnett et al. 
11990] restricted to v re i < 800 km s 1 by a formula 



OCX ( V rel) = AO + £j Gi I" Vre1 )' 

!=1 



(2.16) 



The most r ecent and authoritative compil ation of measurements and calculations was pre- 
sented by \Lindsav and Stebbings B2005I1 . who suggested the following formula for the 
cross section expressed in cm 2 , valid for collision energies E between 0.005 and 600 keV: 



o-cx (E) = 10~ 16 (1 - exp [-67.3/£]f 5 (4.15 - 0.531 ln£) 2 



(2.17) 



A comparison of the cross sections from the 4 formulae is presented in Fig. 12.61 

It is important to note that while all four formulae return similar results for the collision 

speed s relevant for the supersonic solar wind, the recommended formula bv lLindsav and Stebbinss 
1 2005] returns a significantly larger cross sections for low energies relevant for the outer 
heliosheath. Thus adoption of the older formulae may result in a significant underesti- 
mation of the coupling strength between the neutral interstellar gas and the plasma in the 
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2. Solar parameters 
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Figure 2.6: Cross sections for charge exchange reaction between protons and H atoms in 
an en ergy range most impor t ant in the heliospheric physics. The recommended relation 



from \ Lindsav an d Stebbinss 1 2005;! i s com pare d with the formulae used in the past by 



Fite et al. | 1962]; Maher and Tinsley | 1977|] and Warnett et al. Nl990ll 



outer heliosheath, where the secondary population of interstellar H atoms is formed, with 
consequences for the results of heliospheric modeling, as elaborated on in this book in the 
chapter by Izmodenov et al. 



2.3.2 Photoionization 

From solar spectrum to photoionization rate 

Photoionization is a secondary ionization factor of neutral interstellar H, but its sig- 
nificance has recently increased because of the decrease in the average solar wind flux 
observed since the last solar maximum, resulting in a decrease in the intensity of the dom- 
inant charge exchange ionization rate. The rate of photoionization /3 p h (?) at a time t can be 
calculated from the formula 

•'ion 

y6ph (?) = J cr ph (A) F A (A, t) AA (2.1 8) 

o 

where cr p h (A) is the cross section for photoionization for wavelength A and (A, t) is the 
solar spectral flux at the time t and wavelength A; A lon is the wavelength for the threshold 
ionization energy. In the case of hydrogen, the spectral range of the radiation capable of 
knocking out electrons from H atoms is entirely in the EUV range. The cross section for 
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Figure 2.7: Cross s ection for photoionization of hydrogen based on Eq. 12. 191 adapted from 
Verner et al.\ lll996ll . 



photoionization of H can be expressed by the following formula Werner et al. Jl996ll 



o-HphW) = 6.82297 x 10" 



9.3666 4 
VI 



+ 1 



-2.963 



- 2884.69) z /I 



2 -,2.0185 



(2.19) 



where the cross section is expressed in megabarns (Mb), wavelength in nm and A < 
91.18 nm. The cross section is shown in Fig. 12.71 and the importance of various portions 
of the spectrum for the photoionization is illus trated in F ig. [2~8l which shows the integrand 



function from Eq. l2.18l for DOY 122 in 2001 IBochsler et a/.ll201 ll in preparation]. 



It is clear from Fig. I2.8l that the most important for the photoionization of hydrogen is 
the portion of the spectrum immediately next to the photoionization threshold. The pho- 
toionization rate varies because the solar EUV spectral flux varies. Direct measurements 
of the solar EUV spectrum in the entire relevant energy range have been available on a 



2-hourly basis since 2002, when the TIMED/SEE experiment Woods et al.l 1200511 was 



launched. Before, the spectral coverage was intermittent and one had to resort to indirect 
methods. 

Basically, these methods can be grouped into two classes: (1) direct integration using 
Eq. l2.18l with the spectrum Fx (A, t) calculated from proxies, and (2) using correlation for- 
mulae between selected proxies and photoionization rates obtained for the times when the 
spectrum measurements are available. Since measurements covering only a portion of the 
spectrum with a relatively low spectral re solution are also available (SOHO/CELIAS/SEM) 



&Hovestadt et al. , 1995 ; Judse et al. , 1998 ] , a variant of method (2) would require finding a 



correlation between these partial direct measurements and the photoionization rate. None 
of the method gives perfect results; furthermore, a reliable application of method (2) be- 
came possible only after a sufficiently rich database of precise measurements of the solar 
spectrum became available. 

Inspection of Fig. !2.8l shows that the band from 50 to 30 nm, which is the most rele- 
vant for helium and has been observed by CELIAS/SEM, is of secondary significance for 
hydrogen. Since, however, variations in various portions of the solar EUV spectrum are 
correlated to some extent, the evolution of the EUV flux in the CELIAS/SEM bands is a 
reasonable indicator of the evolution of the photoionization rates of H as well. 
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2. Solar parameters 
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Figure 2.8: Integrand function from Eq. l2.18l for the flux obtained from TIMED/SEE for 
2010 DOY 122 and photoionization cross section defined in Eq. l2.19l and shown in Fig. 12.71 



Employment of method (1) requires a proxy model of the solar spectrum in the rele 
vant photon energy range. A number of such models was develo ped in the past, including 
SERF! (also known as HFG orEUV81) iHinteresser et a/.LIl98lll. EUVAC {Richards et all 



1994 |l SOLAR2000 iTobiska et a/.ll200dl . and NRLEUV iWarren et ad[l998al lbl: lLean et 
l2006ll . The models are based on empirical correlations found between var- 
ious portions of the solar spectrum and selected time series of available measurements, 
including typically the F10.7 radio flux, the MgII c / w index, the solar Ca II K index, and 
even the sunspot time series. Methodology and problems of c r eation of s u ch m odels 
have been recently reviewed by \Flovd e t al. [2002, 2003, 200 511: \Lean et al. [2011] and 



'all 



the question of find ing suitable proxies bv \Kretzschmar et al. 1 20061 : XDudok de Wit et al. 



2008, 2009]. Typically, linear correl ations have been so ught for, what does not 



seem to be an optimum solutions at all cases IBochsler et a/.Ll201 



>ugti 
11]. 



Variation in the photoionization rate of H in the ecliptic plane 

The process of ph otoionization of heliospheric hydrogen was extensively discussed 
b\ \Oeawa et al. 1 1995 1 and we will not repeat such a discussion in this Chapter. We will 
briefly present current views on the rate of photoionization of heliospheric hydrogen based 
on actual measurements and models. 

We have developed a model of evolution of the Carrington rotation averages of the hy- 
drogen photoionization rate s that is based on direc tly measured solar spectra from TIMED/SEE 
and a hierarchy of proxies IBochsler et ali l201li in preparation] . From the spectra at full 
time resolution available from 2002 until present we calculated the photoionization rates 
using Eq. 12.181 Since the time series obtained showed clear signatures of flares and local 
particle-background contamination, we filtered it against the outliers beyond two sigmas, 
which show up in the time series of change rates (J3(ti+i) — B{ti)) / (f, + i - f,). From the 
filtered time series we calculated monthly averages, which are shown in the blue color in 
ther lower panel of Fig. 12.91 

Since the coverage by TIMED/SEE is limited in time and intercalibration or compari- 
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son of various measurements in the heliosphere taken at different times requires a knowl- 
edge of homogeneously-derived time series of ionization rates, we used a hierarchy of 
proxies to extend the directly-obtained photoionization rates backwards in time until the 
end of 1947, when the measurements of the F10.7 flux became available. 

We started from directly integrated photoionization rates calculated from the filtered 
TIMED spectra, which now cover the full interval from solar maximum to solar mini- 
mum. We calculated a time series of monthly averages and we found a correlation formula 
between these values and Carrington rotation averaged measurements in Channel 1 and 
Channel 2 from CELIAS/SEM and the time series of Lyman-alpha flux from LASP, which 
is the following: 

/3 Hp h = 5.39758 x 1(T 2() / tot + 2.36415 x 1(T 16 SEM^f 5549 + 

+ 3.98461 x 1(T 16 SEM^f 549 + 2.05152 x 1(T 8 (2.20) 

Using this formula, we calculated the Carrington averages of photoionization rates for 
the entire interval for which the SEM data are available. For the times when they are 
unavailable, but the MgII c / w index from LASP is available, we use another correlation 
formula: 

j8 Hp h = 3.56348 x 10" 6 MgII c/w - 8.5947 x 10~ 7 (2.21) 

When the MgII c / w index is unavailable, we use the following correlation formula with the 
F10.7 flux expressed in the sfu units (i.e. 10~ 22 W nT 2 Hz -1 ): 

/3 Hp h = 1.31864 x 10~ 8 F l > 4 ™ 344 - 1.7745 x 10~ 8 (2.22) 

It is worth pointing out that the exponent at the F10.7 flux is close to 1 /2, not 1, as frequently 
adopted. 

The results of the model with color coding of the sources used is shown in the lower 
panel of Fig. |2.9| for the time interval since 1990 until present. 

The photoionization rate obtained in the way described above can be compa red with 



the ra te from the SOLAR2000 model, used extensively in the previous studies \e.s. \B7.owski et a 



2009]. The comparison is shown in the upper panel of Fig. 12.91 The two models agree to 



about 10 - 15%, with the direct integration model giving almost always higher values. 
Such an accuracy is basically equal to the present accuracy of the EUV measurements, 
especially in the low-energy portion of the spectrum which contributes most to the ioniza- 
tion. 
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2. Solar parameters 



Monthly photoionization rates of hydrogen: 
SOLAR2000 vs present model 




1960 1970 1980 1990 2000 
time |y| 

SOL2000 vs direct integration relative difference 




Monthly photoionization rates of hydrogen: 
direct integration model since 1990 




Figure 2.9: Comparison of the monthly averages o f hydrogen photoionization rate ob- 
tained from the SOLAR2000 uobiska et all l200Ll gray] and direct-integration model 
(blue), exte nded using the SEM/L y-alpha (dark blue), Mg II C / W (green), and F 10.7 prox- 
ies (black) XBochsler et all 1201 ll in preparation] in the upper panel; relative difference 
between the SOLAR 2000 and present model (middle panel), and exploded view of the 
photoionization rate from the present model during the two past solar cycles (lower panel). 
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Latitude variation of photoionization rate 

Similarly as in the case of solar Lyman-alpha radiation, al so radiation in t he pho toion- 
ization spectral region is expected to vary with heliolatitude. \Auchere et al.\ j2005bll con- 
structed a model of equatorial and polar flux in the solar 30.4 nm line, mostly responsible 
for ionization of helium, and demonstrated that a ~ 0.8 pole-to-equator ratio (fluctuating) 
can be expected. A similar modeling for the spectral r ange relevant for photo ionization of 



hydrogen is not available, but coronal observations bv lAuchere et al.\ M2005all suggest that 



such an anisotropy can indeed be expected and that some north-south asymmetry cannot 
be ruled out. 

In lack of a complete model and sufficient data we surmise that a latitude variation 
of hydrogen photoionization rate approximated by a formula similar to Eq. 12.41 may be 
tentatively adopted. The subject certainly needs further studies. 



2.3.3 Electron ionization 



Significance of the electron-impact ionization reactio n for the distribution of neutral 
interstellar gas in the inner heliosphere was pointed out bv lRuciriski and F ahr lll989ll . who 
developed a model of the electron ionization rate based on the local electron temperature 
and density. Further insight into the problem of electron- impact ionization of neutral in 



terstel lar hydrogen inside the heliosphere can be found in lBzowski et al\ [2008]; B zowst 
I 2008h . 

The ionization rate in the electron-impac t reaction at a locatio n described by the radius 



vector r can be calculated from the formula IQwocki et all 1983] 



ySei (r) 



CO 

_ 8tt r 

ml J 



o- e] (E)f e (E,r)EdE, 



(2.23) 



where cr e \ is the energy-dependent reaction cross section, E collision energy and £i on the 
ionization threshold energy, for hydrogen equal to ~ 13.6 eV. Ionization occurs for H atom 
- electron collisions with energies exceeding the limiting energy E loa . Practically, almost 
entire energy of electrons in the solar wind at a few AU from the Sun is in the thermal 
motions. Kinetic energy of an electron moving at a typical solar wind expansion speed of 
440 km s 1 is about 0.5 eV, which is much less than the ionization threshold. Since the 
temperature of the core electron distribution function at 1 AU is ~ 10 s K and of the halo 
component on the order of ~ 10 6 K, the thermal speeds of the solar wind electrons are on 
the order of a few thousand of km s _1 , which is much more than the expansion speed. 
Thus, the expansion speed of the electron fluid can be neglected in the calculation of the 
electron-impact rate. 

The formula for the cross section for electron ionization was proposed bv lLod jl967bll 
and simplified for H by Lotz\ | 1967all : 



ft af g | .ln(£) (l-frexp[-c,-(f -1)]) 

el (E) = 2^ 



EP, 



(2.24) 



where Na is the number of electrons in the ion and the summation goes over the partial 
cross sections for knocking out all individual electrons from the ion. P, is the ionization 
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2. Solar parameters 



Cross section for H electron-impact ionization 
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Figur e 2.10: Cro ss section for electron-impact ionization of hydrogen, defined in Eq. 12.231 
as a function of H atom-electron collision energy in eV. 



potential for a given charge state of the ion, E is the impacting electron energy, a,, , c, are 
parameters specific to a given ion and its charge state and is the statistical weight. For 
hydrogen, there is only one electron to be knocked out and Eq. l2.24l takes the form: 



cr el (E) =4.0 1 -0.60exp 



-0.56 



,13.6 



ln(£/13.6) 
13.6£ 



(2.25) 



where E > 13.6 is expressed in eV. It is claimed by Lotz | 1967a] to be accurate to ~ 10% 
and shown in Fig. 12.101 

Electron density can be adopted from quasi-neutrality and continuity conditions in the 
solar wind and calculated as equal to the local proton density + 2 x local alpha-particle 
concentration: 



= «,(l+2&), 



(2.26) 



where % a is the local alpha-particle abundance relative to solar-wind protons. 

The temperature behavior is more complex. Since electrons are subsonic while pro- 
tons are supersonic on one hand, and since collisions are infrequent on the other hand, 
the electrons cool off faster than protons, and not adiabatically. The distribution function 
of electrons in the solar wind and its evolution with solar distance is fairly complex and 
requires further studies, especially that measurements perform ed on Ulysses using two dif- 
ferent techniques: Q uasi Thermal Noise jlssautier et a/.ll200lll and particl e-measurements 
I Salem et al. , 2001 1 return somewhat discrepant results jlssautier , l 2009ll . Basically, the 
electron distribution function can be decomposed into three components: a warm core, a 
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hot halo (both approximated by the Maxwellian function: 



f e (r, T e ) = n e 



In k B T e (r) 



3/2 



exp 



k B T e (r) 



(2.27) 



where k B is the Boltzmann constant, m e electron mass, E kinetic energy of electron, T e 
temperature of the electron fluid, and r radius-vector of the locat ion in space), and a fluc- 
tuating strahl, stretched along the local magnetic-field direction jPilipp et al. . ll987b[ laFl 
The contribution of the halo population to the net r ate is on the level o f a few percent and is 



an inc reasing functi on of the heliocentri c distance \Maksim ovic et al. . 2005 : Stverdk et al. 



2009]. Estimates bv Wzowski et al.\ 1200811 show that at 1 AU, the ionization rate due to the 



core population of the solar wind electrons is equal to about 0.4 x 10 s and to the halo 
population to less than 0.04 x 10~ 7 s _1 . The amplitude of fluctuations in the electron ion- 
ization rate may reach an order of magnitude, which is much more than the long-time vari- 
ations related to var iations in solar activity. On the other hand, the electron data from Wind 
jSalem et ali l2003ll imply an in-ecliptic solar minimum (1995) rate of ~ 0.68 x 10~ 7 s _1 
and a solar maximum (2000) rate of ~ 0.73 x 10~ 7 s _1 . Thus assuming a constant rate over 

the solar cycle is a reasonable appro ximation. 

Observations done with Ulysses jPhillips et al. , 1995at Issautier et al. , 1998 ; Le Chat et al. 



201 111 suggest that the electron ionization rate is a 3D, time dependent function of the solar 
cycle phase. Both th e tempera t ure magnitude and the cooling rate differ between the fast 
and slow solar wind. \Bz.owski ll2008ll adopted the following radial profiles of the core T c 
and halo T;, temp eratures and the halo-to-core density ratios She = n/,/« f for the slow [after 



Scime et aZ.1, 119941: 



he 



T c = 1.3 -lOV - 85 
T h -- 
0.06 r-°- 25 



.2-10 5 r-°' 38 



(2.28) 



and fast \after \lssautier et ali 1998; Maksim ovic et al. , 2000] solar wind 



T h /T c 

She 



7.5 -lOV ' 64 

13.57 

0.03 



(2.29) 



In both solar wind regimes, the core (n c ) and halo («/,) densities are calculated from the 
equations: 



1 + 2g g 

1 +£hc 
She «c 



(2.30) 



with the solar wind alpha abundance $ a = 0.04, identical in both fast and slow wind 
regimes on average, and n p the proton density. 



'Note that the electron distribution function was recently approximated by Le Chat et al. 1 2010, 201 1] by the 
kappa function, which naturally covers both the core and halo components 
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2. Solar parameters 



Rucinski and Fahr 1 1989] inserted the formulae from equations 12 . 241 l2~26l and !2.27l to 



the integrand function in Eq. l2.23l and obtained a formula for each of the terms i contribut- 
ing to the total electron-impact cross section, which we present here in a slightly modified 
form: 



ySei.i (T e ) = n t 



m e Pi V t ksiT, 



r o 



P . \ r(0,c, + ^)exp[c,] btPi 



k n T f 



Pi + Cjk B T e 



(2.31) 



where T(a,z) is the incomplete gamma function and for hydrogen i — 1, q, ■ - 1, P, ■ - 
13.6 eV, Oj = 4.0 x 10~ 14 , b, = 0.60, a = 0.56 . 



Us ing these relations. LBzows fa [ 2008] employed the approach proposed bv Ltoc;»sfc; am/ Fahr 
I 1989h . assuming tha t the core and halo te mperatures are isotropic (which is not exactly 
the case, as shown b\ \Stverdk et al.\ M2008I1 ). and calculated radial profiles of the electron- 
impact ionization rates separately for the fast and slow solar wind and approximated the 
results by the following phenomenological formulae: 



jSei.f (r, rip) = 



-r exp 



lip 

~2 eX P 



In r (541. 69 In r- 1061.32) + 1584.32 



(In r - 29.17) ((In r - 2.02) In r + 2.91) 
In r (348.73 lnr - 917.39) + 2138.05 



(2.32) 



(2.33) 



(In r - 18.97) ((In r - 2.53) In r + 5.74) 

As is evident in these formulae, the electron-impact ionization rates are parametrized by 
local proton densities n p normalized to 1 AU and are fixed functions of heliocentric dis- 
tance, which differ appreciably from the 1/r 2 profiles that are typical of the solar wind flux 
and photoionization rate, as shown in Fig. 12. Ill The behavior of electron tempe rature in 



the int ermediate solar wind resulting from stream-stream interactions is unknown. Wzowsk 



1 2008] approximated the electron temperature in the solar wind of intermediate velocity 
as a mean value between the fast and slow wind specific temperatures, weighted by the 
percentage of the solar wind density in transition from the slow to fast wind conditions. 

A consequence of the different magnitudes and radial behavior of the solar wind elec- 
tron temperature between the slow and fast solar wind is a pronounced latitude anisotropy 
of the electron ionization rate throughout the solar cycle with an exception for a brief in- 
terval at solar maximum, when solar wind becomes almost spherically symmetric. The 
significance of the electron impact ionization in the overall balance of the contributing 
ionization reactions increases towards the Sun. Eventually, the role of electron ionization 
and its anisotropy is greater in the downwind hemisphere, where neutral interstellar gas 
has already passed the Sun (i.e., the streamlines of the H gas flow have passed their peri- 
helia). Because of the fast cooling of the electron fluid with solar distance, the significance 
of the electron ionization reaction for the distribution of neutral interstellar hydrogen in the 
upwind hemisphere at the distances above ~ 2 AU (i.e., at and beyond the distance of the 
Maximum Emissivity Region (MER) of the heliospheric Lyman-alpha backscatter glow) 
is almost negligible. 



2.4 Evolution of solar wind in the ecliptic plane 

The existence of solar wind was predicted on a theoretical basis bv \Parke~r 1 1958 1 and 
discovered experimentally once the spacecr aft Lunnik II and Mariner 2 left the magneto- 
sphere at the very beginning of the space age IGri neauz et al. , 1960; iNeueebauer and SnvdeA 
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Figure 2.11: Normalized radial profiles r 2 f3 e \ (r) of the equatorial (red) and polar (blue) 
rates of electron-impact ionization of neutral inter stellar H atom s in the supersonic solar 
wind, defined by Eq l2.32l and ( 12.33b [adapted from \BzowskL 1200811 respectively. 



Regular measurements of its parameters begun in the beginning of 1960s and up to 



now, data from more than 20 spacecraft are available, obtained using various techniques 
of observations and data processing. 

The highly supersonic solar wind in the ecliptic plane consists in fact of a sequence of 
various interleaved components: a "genuine" slow solar wind, fast solar wind, solar wind 
plasma from stream-stream interaction regions, and (intermittently) interplanetary coronal 
mass ejections (CME). The Mach number of the flow varies from ~ 3 to ~ 10. The balance 
between the populations changes with solar activity. 

Solar wind expands almost radially (i.e., its velocity has very small non-radial compo- 
nents in comparison with the radial component), with the speed basically invariable with 
solar distance between the outer boundary of the acceleration region near the Sun (located 
inside a few solar radii ) and approximately 10 AU, where the effects of mass loading due 
to the ionization of neutral interstellar gas \Fahr and Rucinski, 1999; Fahr and Rucinski, 



2001, 2002] become measurable \Richardson et al., 1995, 2008]. The overall slowdown 



continues up to a few AU upstream from the termination shock, where the flow speed, 
by then already reduced by about 60 km s _1 relative to the speed at 1 AU, i s additionally 



slowed down by the componen t of protons reflected at the termination shock jLiewer et al. 



1993; Richardson et al., 2008]. For the purposes of modeling of the neutral gas distribution 



observed from 1 AU the reduction in solar wind speed with distance may be neglected. 
Outside the acceleration region of a few solar radii, solar wind density falls off with 
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Figure 2.12: Ionization rates in the ecliptic plane from all relevant processes (black). 
Blue marks th e charge exchange rate (calcu lated from Eq. 12.141 using the cross section 
from Eq. 12. 171 \ Lindsay and Stebbingsl 120051). green the photoionization rate, orange the 
electron-impact ionization rate. 



solar distance as 1/r 2 and this relation, stemming directly from the continuity equation, is 
not significantly altered by the interaction of the solar wind plasma with neutral interstellar 
gas. Indeed, with the main reaction being charge exchange, a reaction act results in a shift 
of the proton in phase space (a core solar wind proton is replaced with a pickup proton 
moving at a different velocity), but not in a change in the local density. Only the two 
secondary reactions, photoionization and ionization by electron impact, actually inject new 
protons into the pickup ion region in phase space, thus increasing the local concentration 
of protons. Since the two secondary ionization processes contribute only ~ 20% of the 
total ionization rate at 1 AU, as it will be clear later in this chapter, and since the PUI 
content in the solar wind at the termination shock is on a level about 20%, the net increase 
in the absolute density of solar wind above the level dictated by the 1 /r 2 relation can be 
estimated as an excess of 20% x 20% = 4% , which is negligible in comparison with 
the uncertainty in the solar wind density, which will be discussed in a further part of this 
chapter. 

The probability distribution function (PDF) of hourly averages of densities is log- 
normal (i.e., logarithms of densities are distributed normally), while velovity features a 
complex PDF which can be approximately described as a sum of at least two normal distri- 
butions with differently l ocated peaks of different widths. They are illustrated in Fig. 12.131 
[see also ltfe Toma , 1201 I I. 

From the viewpoint of charge exchange rate, details of the PDF of solar wind pa- 
rameters are not as important as their mean values at appropriate time scales. Since the 
parameters of solar wi nd vary with the solar cycle and, as it was recently discovered by 
McComas et al. j2008tl . also secularly, and since neutral interstellar hydrogen, which is 



the background for the heliospheric Lyman-alpha backscatter glow, features variations 
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Figure 2. 13: Histograms of hourly averages of solar wind density (adjusted to 1 AU) (left) 
and speed (right) from the entire OMNI-2 database, which currently spans an interval from 
1963 until 2011. 



on time scales comparable to the solar cycle length I Rucinski and Bz.owski , 1995bll . long 
and well-calibrated series of solar wind parameter measurements is a very important asset 
for heliospheric studies. Historically, measurements of solar wind speed obtained from 
various experiments generally agreed among themselves with an accuracy of ~ 5%, but 
systematic differences between density values from different experiments up to ~ 30% ex- 
isted. Hence any study of a long-term behavior of solar wind required intercalibration of 
the results from different experiments. Such an initiative brought the OMNI data collec- 
tion fKing and Papitashvili a vailable at at http : //omniweb . gsf c . nasa . gov/| [see also 
\Kine and Papitashvilh\2Q05\ . where historical and present measurements of the solar wind 
density, velocity, temperature, alpha abundance, and magnetic field vector were brought to 
a common calibration. 

Originally, the OMNI collection was created in the 1970s by the National Space Sci- 
ence Data Center (NSSDC) at the NASA Goddard Space Flight Center. In 2003, a succes- 
sor database OMNI-2 was made available, which has been maintained until present. The 
OMNI-2 data collection was basis for the present study of solar wind parameters in the 
ecliptic and is used in this chapter to construct the charge exchange ionization rates. 

The ea rly period of t he OMNI databases ( 1963 to 1971) includes data from multiple 
spacecraft \Bonetti et al\ II 196911 : Weusebauen 1197011 . the middle-period data are mostly 
from IMP-8 and span from 1971 to 1994, a nd the later per iods, from 1994 until p resent, in- 
clude mostly data from IMP-8, Wind/SWE jKaspeAWO^ . and ACE/SWEPAM jMcComas et al. 



1998]. Since there was no overlap between the early and middle period, the data from the 



early period in the OMNI-2 database are adopted unchanged from the original OMNI col- 
lection. The data from the early period were extensively intercalibrated between various 
spacecraft doing observations during that time, but no intercalibration with the middle and 
recent periods was possible. Still, owing to the overlap of the data between the middle and 
late periods, it was possible to perform an intercalibration between the data from these two 
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periods. 

Most of the solar wind plasma data used in the OMNI collection were obtained from 



the MIT Faraday Cups [Bridge et al, 1965; Lyon et al., 1968; Lyon et al., 1967; Lazarus and Paularena, 
1998 1 and LANL electrosta t ic analyzers iBame et al. . \91 \ . Hundhausen et all 19671: Osilvie et al . 



1968; Feldmanetal., 1973; \Asbridee et alW\91(k\Bame et al\. \\91?,A U \McComas et al.l 
1998I1 . In the OMNI-2 series. \Kins and Papitashv ili [2005] adopted measurements of so- 
lar wind parameters performed by the Wind spacecraft as basis for the common calibra- 



tion. T hey followed in this respect an analysis performed by Kasuer [2002] ; Kasuer et al. 



1 2006]. The preparation of the data published in the OMNI collection involves removing 
from the original high-resolution data supplied by the Principal Investigators of the exper- 
iments possible Earth bow shock contamination and incomplete records, and subsequently 
time-shifting the data from the spacecraft location to the Earth location. The primary 
source of the solar wind data are nowadays Wind measurements, which are ultimately 
replaced with measurements from ACE, which unlike the data from Wind, are free from 
possible bow shock contamination because the ACE spacecraft operates near the Earth at a 
Lissajous orbit close to the LI Lagrange point about a million kilometers upstream of the 
Earth bow shock. Because of different time scales of the processing of the data from vari- 
ous spacecraft, typically an interim data product becomes available once the first data are 
obtained, which are superseded with the final product when all the data needed become 
available or are declared as unavailable. This results in some changes in the published 
time series during the time. Our experience shows, however, that the changes are seldom 
significant. 

The data from different experiments are scaled to a common calibration using linear 
fits based on results of linear regression analysis. It is also important to mention the signif- 
icance of the correlation work that the OMNI team performed on the data from the IMP-8 
spacecraft and early Wind and ACE measurements. IMP-8 operated for 28 years in 1973— 
2001 interval and provides a bridge between the early and present observations, enabling 
presentation of a more or less homogeneous series of solar wind parameters shown further 
down in this chapter. 

The result of the intercalibration process is a time series of hourly-averaged solar wind 
parameters. Because of the varying quality of individual records, the time coverage of the 
parameters is not uniform and gaps may exist in some parameters, while correct data for 
the same time interval may be available for other. Since the distribution function of solar 
wind is inhomogeneous and varies rapidly in time, the values of solar wind parameters 
retrieved from observations depend on the method used to process the data. Typically, the 
LANL team take moments of the observed distribution function to calculate the density, 
speed and temperature of the solar wind, while the MIT team fit the measurements to an 
anisotropic Maxwellian or bi-Maxwellian function using a nonlinear fit method. To assess 
differences resulting from the two aforementioned approaches, the MIT team calculated 
the density, speed and temperature from the Wind/S WE distributions using both mentioned 
methods. 

King and Papitashvili at http://omniweb.gsfc.nasa.gov/ extensively discuss the differ- 
ences and correlations between data from various sources. They show that the velocities 
are very well linearly correlated, with the coefficients of the relation vswe - a + b vace 
equal too = -2.135 + 0.387 , b - 1.010+0.001. In the case of densities, it is the logarithms 
of density which are linearly correlated and the coefficients of the formula log «swe = 
a + b log«ACE slowly vary with speed, a changing from 0.006 for v < 350 km s _1 to 
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0.091 for 350 < v < 450 km s _1 to 0.082 for v > 450 km s _1 and b changing from 1 for 
v < 350 km s _I to 1.036 for 350 < v < 450 km s _1 to 1 for v > 450 kms 4 . This brings 
up differences on the order up to 20% between ACE and Wind, which are comparable to 
the uncertainty in density coming up from the application of various methods of parameter 
derivation discussed earlier. In a nutshell, while a very good correlation of speeds is ob- 
tained, the correlation between logarithms of densities is close to linear, but with a scatter 
of approximately 30% around the fit line. This is probably a good measure of inherent 
uncertainty of the densities even without uncertainties in the absolute calibrations. 

Based on the OMNI database, we constructed a time series of Carrington period- 
averaged parameters of the solar wind normalized to 1 AU, with the grid points set pre- 
cisely at halves of the Carrington rotation intervals. Small deviations of the times from 
the halves of the rotation periods were linearly interpolated. Averaging over the Carring- 
ton rotation enables constructing of a consistent model of the ionization rate, where axial 
symmetry of the solar parameters is safeguarded. Such a procedure seems to be the only 
reasonable one at present. Solar wind parameters show considerable variations during one 
solar rotation period, with quasi-periodic changes from slow to fast solar wind speeds and 
related changes in density (see Fig. 12. 15b . The time scale of the changes of the fast wind 
streams is comparable to the solar rotation period and thus constructing of a full and accu- 
rate model of the solar wind variation as a function of time and heliolongitude is presently 
not feasible because of the lack of sufficient data. Thus, we developed an equally-spaced 
time series of the solar wind in-ecliptic densities and speeds, which are well suited for 
the purpose of global heliospheric modeling and to calculate the charge exchange rates of 
heliospheric neutral atoms. The time series of density, speed, and charge exchange rate in 
the approximation of neutral H atom stationary relative to the Sun is presented in Fig. l2.14| 

The time interval shown in this figure starts before the solar activity minimum in 1986 
and includes the solar minima of 1995 and 2007, as well as the two recent maxima of 
1990 and 2001. One observes a striking difference in the appearance of the solar wind 
equatorial parameters in comparison with the behavior of solar EUV radiation. Neither 
density nor speed seem to be correlated with the level of solar activity. There is no clear 
minimum-maximum-minimum variation, which was clearly seen in the EU V-related time 
series. Speed shows multi-timescale variations, but its mean value is basically unchanged 
in time. By contrast, density features a secular change, which begun just before the last 
solar maximum and leveled off shortly before the present minimum. The overall drop in 
the average solar wind density is on the order of 30% between 1998 and 2005; it seems 
that the drop in the charge exchange rate is even more dramatic, namely by a factor of 
two. The present rate of charge exchange is at a level similar to th e charge exchange rate 



observed by Ulysses a t the poles during its first fast latitude scan tPhillips et all 11995a: 



McComas et al. . Il999h . Thus the solar wind density features a "plateau" until 1998, then 



a "cliff" and a "foot" starting from 2002. Within the "foot" there seem to exist fluctua- 
tions of density and anticorrelated fluctuations of speed. These can be associated with the 
persist ence of coronal holes at equatorial latitudes, as convincingly illustrated bv \de Toma 
|2011]. This overall long-standing drop in the solar wind flux results in an overall en- 
hancement in the neutral interstellar gas density near the Sun. The secular variations in 
the solar wind, on time scales significantly longer than the solar cycle time scale, suggest 
that the heliosphere does not evolve strictly periodically and that monitoring of these vari- 
ations is an essential element of any effort aimed at quantitative analysis of all kinds of 
heliospheric observations. 
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Figure 2.14: Carrington rotation averages of solar wind density adjusted to 1 AU (upper 
panel), speed (middle panel), and flux (lower panel) calculated from hourly averages from 
the OMNI-2 database jKine and Papitashvili 1200511 . The thin vertical lines mark the time 
interval shown at daily resolution in Fig. 12.151 
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Figure 2.15: Daily averages of solar wind density (adjusted to 1 AU) (upper panel) and 
speed (lower panel) calculate d from hourly averages obtained from the OMNI-2 database 
I Kins and Papitashvill 1200511 for an example time interval covering 7 full Carrington rota- 
tion periods, presented to illustrate the complex structure of solar wind parameter evolution 
at different time scales and the approximate anticorrelation of density and speed. 
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2.5 Latitudinal structure of the solar wind and its evolu- 
tion during the solar cycle 

Shortly after the discovery of solar wind, a question whether or not it is spherically 
symmetric appeared. It is currently known to feature a latitudinal structure which varies 
with the solar activity cycle. While direct observations of the solar wind in the ecliptic 
plane have been conducted since the early 1960s, information on its latitudinal struc- 
ture was avai l able only fr om indirect and mixed observations of the cometary ion tails 



^Brandt et al. . 1972 , 1975 1. The situation changed when radio-astronomy observations of 
interplanetary scintillation and appropriate interpretation of observations of the Lyman- 
alpha helioglow became available. Apart from the in situ measurements obtained from 
Ulysses, these two techniques remain the only source of global, time-resolved information 
on the solar wind structure. 



The launch of the Ulysses spacecraft IWenzel et al. . 1989 1. the first and up to now the 



only interplanetary probe to leave the ecliptic plane and sample interplanetary space in 
polar regions, improved our understanding of the 3D behavior of solar wind by offering 
a very high resolution in latitude, but a poor resolution in time (the same latitudes were 
visited only a few times during the ~ 20-year mission). Hence the studies of solar wind 
parameters as a function of time and heliolatitude are still a work in progress and therefore 
are discussed in a separate section in this chapter. 



2.5.1 Historical perspective: insight from interplanetary scintillation 



Interplanetary scintillation measurements involve radiotelescope observations of re- 
mote compact radio sources (like quasars) and searching for fluctuations of the signal that 
come up because of diffraction of radio waves on electron density fluctuations ~ 200 km 
in diameter occurring in the line of sight. Specifically, the observable is the scintilla- 
tion index, defined as a quotie nt of the r.m .s. of the observed intensity fluctuation to 

the mean intensity of a source \ Manoharan[ 1993 1. Owing to a correlation between the 
fluctuations amplitude and solar wind speed it is possible to deduce the solar wind speed 
fro m careful analysis of registered diffraction patterns in the observed radio source sig- 



nal iHewish et ali 1196411 . The scintillation signal is the sum of waves scattered along 
the line of sight to the observed radio source. Most of the scattering occurs at the clos- 
est distances to the Sun along the sightline because the strength of the electron density 
fluctuations rapidly decreases with solar distance. Had the solar wind been spherically 
symmetric, it would be possible to define a weighting function and the IPS "midpoint" 
speed would be the spatial average of th e solar wind speed centered at the closest point of 
the sightline I Coles and Maasoe . 1972 1. But solar wind is not spherically symmetric and 
features streams, which result in discrepancies of the measured speed from the actual one 
II Houm inen 1 1 97 ill . 

An improvement in the accuracy of measurements of solar wind speed can be achieved 
owing to a better identification of the scintillation patterns in the sightlines when they are 
simultaneously observed by multiple stations displaced longitudinally on the Earth surface , 



and by using the solar wind tomography te chnique Uackson et all Il997t Uackson et al. 
\\99mKonma et alW\99&\Asai et a/.LIl998ll . 

Despite the high sophistication of the tomography technique, its accuracy inevitably 
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depends on purely geometrical considerations (the telescopes are located in the north hemi- 
sphere of the Earth, whose orbit is tilted at an angle of 7.25 deg to the solar equator), on the 
number of observations available, and ultimately on the fidelity of the adopted correlation 
formula between the scintillation index, which is directly measured, and the solar wind 
speed. W hile correlating these quantities was possible early enough for the equatorial so- 
lar wind lHarmon. \l975]. the out-of-ecliptic IPS measuremen t s coul d only be calibrated 
once in situ data from Ulysses became available \Koiima et al, 2001]. Even then one has 
to keep in mind that ideally, such a calibration should be repeated separately for each solar 
cycle because, as discussed earlier in this chapter, solar wind features secular changes. 

Early meas urements of solar wind speed using the IPS technique brought mixed con- 
clusions: while Dennison and Hewish 11967] di scover ed an increase in solar wind speed 
outside the ecliptic plane. IHewish and Svmonds 111 96911 did not fi nd such an increase. Fur- 
ther observations, however, reported \sv \Coles andRicketl lll976ll . clearly showed that the 
solar wind is structured, with a band of slow speed around the solar equator and polar caps 
of a much faster wind. 

An extensive program of IPS observations of solar wind, initiated in 1980s in the Solar- 
Terre strial Environmental Laboratory at the Nagoya University (Japan) \Konma and Kakinuma , 
1990], resulted in a homogeneous dataset that spans almost three solar cycles and en- 
ables studies of the evolution of solar wind speed profile with changes in solar activity 
\ Konma and Kakinuma. 1 1 9871 . Even before the introduction of the Computer Assisted 
Tomography technique, they suggested that the solar wind structure varies with solar ac- 
tivity, with the slow wind reaching polar regions when the acti vity is high. Supported and 



interpreted by the tomography technique mavashi et al. . l2003h . it enables detailed studies 

activity conditions \Koiima et q/.Ul999t 



Koiima et al. 


2001 


Ohmi et all\lM\\ 2003 


Koiima et al. 


2007 


; Tokumaru et al.. 2009. 



Fuiiki et ali l2003bl 0, 0; \Ohmi et al. ., 2003; 
2010]. The solar wind speed data obtained 



from these observations will be discussed later in this chapter. 



2.5.2 Historical perspective: insight from heliospheric backscatter 
glow 

Observations of the Lyman-alpha heliospher ic backscatter glow, carried out since the 
beginning of 1970s IBertaux and Blamoni 1 1 97 ill . have also been used as a tool to discover 
the 3D structure of solar wind and its evolution with solar cycle. The bi-modal structure 
of the solar wind, with the slow and dense wind in an equatorial band and rarefied, fast 
wind at the polar caps, results in a distinctly structured charge exchange ionization rate. 
Since the direction of inflow of neutral interstellar gas is very close to solar equator, the 
inflowing atoms whose orbits happen to be in a plane close to solar equator spend their 
entire travel through the heliosphere in the region of increased probability of ionization. 
In contrast, those traveling in planes inclined at greater angles to solar equator spend rel- 
atively little time at low heliolatitudes or even do not get close to the region of increased 
ionization at all and hence more of them are able to survive the travel towards the Sun 



[see Fig. 1 in Xhallement et all 1985al for an illustrative sketch]. As a result, a region of 
reduced density of neutral interstellar hydrogen gas is created close to solar equator. At 
heliolatitudes farther from the equator the density of the gas at similar distances is higher. 
The gas is illuminated by an approximately spherically symmetric solar Lyman-alpha ra- 
diation, which is backscattered by resonance fluorescence. Since at equatorial latitudes 
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the density of the gas is reduced and the illuminating flux is almost homogeneous in he- 
liolatitude, the intensity of the backscattered radiation is lower at equatorial latitudes than 
at the polar caps. This equatorial dimming of the helio glow, referred to as the heliospheric 



groove, was observed already in 1970s and 1980s | Kumar and Broadfook 1 19781 1 19791: 



Lallement et al. , 1985bi 1986 : Lallement and Stewart , 1990ll and correctly interpreted 



as 



due to the enhanced ionization level at equatorial latitudes because of the anisotropy of 
solar wind. 

It is important to stress, however, that analysis of the heliospheric Lyman-alpha backseat - 
ter glow is only able to yield the latitude structure of the total ionization rate of neutral 
interstellar hydrogen. From the view point of density structure of neutral interstellar hy- 
drogen near the Sun the nature of the ionization processes is not important, only the result 
they produce, i.e. a decrease in the total density at an equatorial latitude band. Hence no 
differentiation between the charge exchange, photoionization, and electron-impact ioniza- 
tion can be made based solely on the heliospheric glow analysis. Since, however, charge 
exchange is the dominant process and since photoionization is only slightly anisotropic in 
latitude, in the first approximation the latitude variation in the total rate can be regarded as 
latitude variation in the charge exchange rate, which is proportional to the latitude varia- 
tion in the total solar wind flux, modulated by the dependence of the charge exchange rate 
on s olar wind speed ( see E qs ( 12. 141 and ( 12.171 )). 

| Lallement et al. , 1985all proposed a model of the charge exchange rate between the 



solar wind protons and neutral interstellar H atoms as a function of heliolatitude <p defined 
by the formula 

y8cx(0)=ySo(l-Asin 2 0) (2.34) 

This is a two-parameter relation, where Bq corresponds to the rate at the equator and A is 
a pole-to-equator amplitude, which can be fitted from observations of the helioglow. This 
formula reproduced approximately the limited observations of the groove obtained at the 
early stage of the research. It is able to reproduce various pole-to-equator contrasts in the 
ionization rate, also the situation when the contrast virtually disappears and solar wind 
becomes approximately spherically symmetric. However, the profile of the ionization rate 
obtained from this formula has a full width at half maximum ~ 45 deg and is perfectly 
symmetric about solar equator. Consequently, is able to reproduce neither the north-south 
asymmetries in the solar wind nor the situations when the range of the slow solar wind 
significantly differs from ~ 45 deg. 

Despite these deficiencies, Eq. ( 12.341 ) was successfully used to infer qualitatively the 
solar wind structu re by a number of authors iLallement et ali Il985al \Prvor et ~al\. Il99£ ' 



Prvor et al. , 20031 e.g.]. The conclusions were similar to those obtained from the IPS 



analysis: during solar minimum solar wind is latitudinally structured, with a band of en- 
hanced flux at the equator and two polar caps of a rarefied and fast wind. During solar 
maximum the ordered structure changes and the polar caps become almost fully covered 
with the slow wind. However it became clear that the simple model given by Eq. l2.34l is 
not fully adequate to describe the reality and a need for more observations became evident. 
Bertaux et all ill 995 ] proposed an experiment to study solar wind anisotropies using the 
technique of analysis of the heliospheric Lyman-alpha backscatter glow, which was im- 
plemented in the French/Finnish project SWAN onboard the Solar and Heliospheric Ob- 
servatory mission (SOHO). Already shortly after launch it became evident that the early 
conclusions on the evolution of the heliospheric groove, and thus the solar wind, with the 
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solar activity cycle are supported IBertaux et al. . 1 1 99611 1 9971 1 1 99<A\Kvrold eTali 1 1 998h . 
but the formula used to describe the latitude profile of the ionization rate needs modifica- 
tion. Thus in the latter work, Eq. 12.341 was modified to describe separately the northern 
and southern hemispheres: 

Pcx (<f>) = Po [© (0) (l - A N sin 2 <p) + © {-cp) (l - A s sin 2 <f>) + B (</>)] , (2.35) 

where © is the Heavyside step function, An,A$ are the separate anisotropy parameters for 
the northern and southern hemispheres, and B (<p) is used to narrow the width of the equa- 
to rial band of enhanc ed solar wind flux. Even more sophisticated approach was proposed 
b\ \Summan en 1 1996], who suggested to approximate the equatorial band in the latitudinal 
profile of the total ionization rate by 



for heliolatitudes -40 deg < 
±40 deg, and by 



Pcx (0) = Po [l - A (t) sin 2 (c 0)] (2.36) 
< 40 deg, where c = 9/4 limits the equatorial band to 



Pcx (4>) = Po exp 



0.2P 



(t-P/2) 



(2.37) 



where P is the solar cycle length and t time. In this formula, there was no north-south 
asymmetry allowed, but it was possible to homogeneously reproduce variations of the 
anisotropy parameters in tact with the solar activity cycle. 

Since the north-south anisotropics in the ionization rate and their evolution with solar 
activity were evident on one h and, and on the other hand the first fast latitude scan by 



Ulysses \Plulliv s et q/.L I1995al suggested that the profile of solar wind parameters can be 
approximated by an equatorial plateau (with a "rough surface" of the gu s ty slow wind) 



standi n g out from a fl at "foot" of the fast polar wind, \Bz.owski I 2001all : Bz.owski et al. 
1 20021 : \Bzowskn 1200311 suggested to approximate it by the formula 
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where <p is heliographic latitude and is a shape factor; y6cx,poi is the average ionization 
rate at the poles and the term (/?cx,poi + (>cx <p) describes the north-south asymmetry of the 
polar ionization rates; the term (y8sw, P oi + Sex <f>) + (ySsw.eqtr (0 - jScx.poi) for <p = corre- 

-ln2( V A, : fc ) A 



sponds to the charge exchange rate at solar equator; and the term exp 

describes the latitudinal dependence of the ionization rate [see also \Bz.owski 1200811 . The 
shape of the central bulge is controlled by the exponent ; for N = 2 the shape is Gaussian 
for = 8 it is close to rectangular. 



The parameters in the model bv \Bzowski et al. 1 2003 1 are the north and south bound- 
aries of the equatorial slow wind band 4>n and (p$ , the equator/north pole and equator/south 
pole ratios of the ionization rate, and the polar north/south asymmetry parameters. In order 
to obtain absolute values of th e ionization rate, an in dependent assessment of the i onization 



rate a t the equator is needed. Wzowski et al. 12003] used the theory developed bv lBzowsk, 
1 2003] and interpreted a carefully selected subset of SWAN observations. To maintain 
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as much symmetry in the data as possible, and to include in the data the full span of he- 
liolatitudes with a well balanced presence of measurements from all heliolatitudes, they 
chose observations taken within a week of the passage of SWAN through the projection of 
inflow axis ion the ecliptic plane, i.e. at the beginning of June and December of each year, 
and they restricted the field of view to a narrow band going through the pr ojections on sky 



of sol ar equator and poles. They eliminated "searchlights" [discovered bv \Bertaux et al. 



2000], i.e. reflections on sky of the point-like active regions on the solar disk, traveling 
across the sky with the angular velocity of solar rotations, and cleaned the data from the 
contamination by extraheliospheric "chaff" (Milky Way, stars etc.). To eliminate possible 
bias from an imperfect absolute calibration, they tried to fit the lightcurves normalized to 
equatorial values instead of attempting to fit the absolute in tensity of the helioglow. 

In agreement with other studies, wzowski et al.\ 1200311 found that the ordered struc- 
ture of the solar wind present during solar maximum disappears with the increasing solar 
activity. The boundaries between the fast and slow wind regions move polewards and ul- 
timately at solar maximum the slow wind encompasses the whole space. The motion of 
the fast/slow wind boundaries in th e north and sou t h hem ispheres were found to be shifted 
in phase by approximately a year. Wzowski et al] B2008I1 discovered that the areas of the 
po lar fast wind region s are l inearly correlated with the areas of the polar holes observed 
b\ \Harve\ andReceM B2002I1 . which enabled them to calculate the variation in the bound- 
aries between the fast and slow solar wind <p$, <ps for the tim e span of the polar coronal 
holes observations, i.e. from 1990 until 2002. Wzowski et al\ 1200811 suggested also that 
the evolution of (p^, <t>s can be approximated by a formula 



(t) = <f> + 0i exp [- cos 3 (w f)] , 



(2.39) 



where a> p was obtained as 2^/Lyman-alpha main period. The validity range of this ap- 
proximation is limited to the time interval from 1990 to 2002. 

Based on this model of ionization, it was possible to infer the evolution of solar wind 
speed and density as a function of time and latitude. This topic will be covered in a further 
part of this chapter. 



The model from Wzowski et al. [2003]; Wzowski et al. .] ll2008ll was able to reproduce 
more correctly the latitudinal span of the slow wind region and its evolution with solar ac- 
tivity, but was inadequate to correctly reproduce the boundaries themselves. Observations 
from SWAN showed that especially during the transition phases of solar activity, the pho- 
tometric latitudinal profiles of the groove are complex and variable in time . To address this 
proble m, it turned out to be necessary to use an approach proposed by \Summanen et al\ 
[1993] and to model the ionization rate as a multi-step function, with different levels in 
fixed, though arbitrarily sel ected latitudinal bands. Such a model of the ionization rate 
was implemented already b\ \Lallement et a/.llll985all . but then the values in the latitudinal 
"slots" were filled with the values o btained from the analyt i cal models discu s sed earlier 
in this chapter. Now, as discussed bv lOuemerais et al. [2006]; \Lallement et al.\ j2010ll , the 
ionization rates in the latitudinal bands are free parameters fitted to the maps of helio- 
spheric backscatter glow, without any assumptio ns on relations between the neighboring 



bands In contrast to the app roach exercised by Wzowski et al.\ M2003H . \Ouemerais et al. 



I 2006ll : \Lallement et al.\ 11201011 used all the data available, cleaned only with use of a mask 
intended to cut off the known extraheliospheric "chaff". Thus a time series of ionization 
rate profiles at a resolution of ~ lOdeg in heliolatitude and ~ 2 days in time was obtained 
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from fitting of the model to the filtered full-sky maps and subsequent scaling to the in- 
ecliptic ionization rates o btained from in situ m easurements. An illustration of results of 
this analysis is shown i mLallement et al\ 1201011 . 



2.5.3 Historical perspective: Ulysses measurements 

Ulysses, launched in October 1990, has been the first spacecraft to traverse polar re- 
gions of the helio sphere and provide a unique view of solar wind jWenzel et artll989t 



Smith et a/.LI1991ll . After a cruise to Jupiter carried out close to the ecliptic plane, it per- 



formed a Jupiter gravity assist maneuver which cast it away from the ecliptic plane on an 
elliptical nearly polar orbit with an aphelion of ~ 5.5 AU and perihelion ~ 1.4 AU in a 
plane nearly perpendicular both to the ecliptic plane and solar equator as well as to the 
inflow direction of interstellar gas. The period of the orbit was about 6 years. The he- 
liolatitude track of Ulysses is shown in Fig. 12.161 The spacecraft was launched at solar 
maximum and the radial cruise was completed just when solar activity was beginning to 
decrease. The first dive towards the south solar pole was carried out during an interval of 
decreasing activity, followed by the first so-called fast latitude scan, when the spacecraft 
coasted from south to north solar pole, covering almost full span of heliolatitudes during 
just about a year at solar minimum activity conditions. This fast latitude scan preceded 
by about a year the first observations of the helioglow by SWAN/SOHO, which provided 
an opportunity to calibrate the model of ionization rate obtained from analysis of the he- 
liospheric Lyman-alpha glow. Ulysses continued on its polar orbit, performing its first 
full slow latitude scan during an interval of increasing solar activity. The following sec- 
ond fast latitude scan occurred in a totally different solar cycle phase, namely during solar 
maximum. Again this scan took about a year and was performed during an interval of 
dynamically variable solar wind structure. The following slow latitude scan was, in con- 
trast to the first slow scan, during decreasing solar activity and afterwards the spacecraft 
performed its last fast latitude scan, again during a solar minimum activity interval. 

The geometry and timing of Ulysses trajectory brought a unique dataset of direct in situ 



meas urements of pa rameters of solar wind p lasma, obtained from SWOOPS \ Bame et al.l 
1992 1 and SWICS 1 Gloeckler et ali 1992tl experiments. The discoveries and findings 



from the plasma mea s urements were presented in dozens of scientific papers [see, e .g., 

Phillips et ai[\\995\k \k \Marsden and SmithUl99li\McComas et aTlll998i[l^l2000b[ laL l2002al lbl l2003ll2006ll2008h 



An additional benefit from this unique mission is the use of the Ul ysses flight spare plasm a 
instrument with only minor modifications on the ACE mission iMcComas et q/.l Tl9981. 
which facilitates intercalibrating the Ulysses measurements with the OMNI time series. 

The Ulysses solar wind dataset is unique and invaluable because it is the first and only 
direct, in situ measurement of the solar wind parameters outside the ecliptic plane. The 
evolution of solar wind speed, adjusted density and flux during the previous and current 
solar minima and during the previous maximum are compiled in Fig. 12.171 where the pa- 
rameter values are averged over 10-degree bins in heliolatitude. Adjustment throughout 
this text means scaling to 1 AU assuming an average dropoff with heliocentric distance 
as 1/r 2 . It can be seen that the heliolatitude structure during the two minima is basically 
similar, featuring an equatorial enhancement in density with the associated reduction in 
velocity (the slow wind region), and that during solar maximum the slow wind expands 
to all heliolatitudes (see also Fig. 12.201 and discussion of time scales in the variation of 
the solar wind structure). However, the region of slow wind seems to reach farther in 
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Ulysses and Earth heliolatitude vs solar activity 
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Figure 2. 16: Illustration of heliolatitude of Ulysses (red) and Earth (green) during the time 
span of Ulysses mission. The pale blue line is the F10.7 solar radio flux, superimposed to 
correlate variations in solar activity with Ulysses heliolatitude during its more than 3 orbits 
in a polar plane almost perpendicular to the inflow direction of neutral interstellar gas. 



heliolatitude during the last solar minimum than during the minimum of 1995, which, in- 
terestingly, is much less conspicuous in density . A striking feature i s a strong reduction in 
flux, reported already bv \McComas et ali t200^A : \Ebert et au j2009ll . The reduction in flux 
is visible as a continuous trend from the 1995 minimum through the 2002 maximum until 
present. Another int eresting trend is a variati on of ~ 1 km s _I deg" 1 in the fast polar so lar 
wind, discovered bv \McComas et al. 1 2000bll and expanded upon bv \Ebert et all J2009]. 



2.5.4 Retrieval of solar wind evolution: introduction 

In this section we present an attempt to retrieve the history of evolution of solar wind 
density and speed during the past solar cycle to provide the modelers working within the 
FONDUE group with a tool to enable them to develop a model of distribution of neutral 
interstellar hydrogen [see Izmodenov et al., this volume] that could be used as background 
for inter-calibration of various sets of UV observations. We attempt to use all relevant 
datasets, paying special attention to absolute calibrations and possible biasing. As it will 
be shown in the further part of this section, this is still the work in progress. Since an 
essential part of information must be drawn from a careful interpretation of the Lyman- 
alpha helioglow observations from SWAN/SOHO, which are subject to improved absolute 
calibration (which is one of the goals of the FONDUE group), the effort of constructing a 
homogeneous set of solar wind parameters must be an iterative one. 

In the following, we will present a construction of the procedure to retrieve the vari- 
ability of solar wind speed and density in time and heliographic latitude with use of the 
available datasets. The procedure to retrieve the time and heliolatitude evolution of solar 
wind relies on the absolute calibration of the OMNI dataset both in speed and density in the 
ecliptic plane. Out of ecliptic, the baseline is the absolute calibration of Ulysses measure- 
ments and interplanetary scintillation observations (IPS). The prime source of information 
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Figure 2.17: Solar wind speed (upper panel), adjusted density (middle panel), and ad- 
justed flux (lower panel) as a function of heliolatitude for the first (blue), second (red) 
and third (green) Ulysses fast latitude scans. The parameters are averaged over 10-degree 
heliolatitude bins. 
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Sun 

Figure 2.18: Illustration of geometry of the line of sight (LOS) in an attempt of remote 
sensing observation of a solar polar region. The observer is close to the Earth in the 
ecliptic plane and aims its instrument (e.g., a radio telescope antenna or a Lyman-alpha 
photometer) at a target so that the line of sight crosses a cone with a small opening angle 
centered at the pole. The signal is collected from the full length of the line of sight, but the 
contributions from various parts are different and depend on the observations technique. 
In the case of IPS observations, the strongest contribution to the signal is from the point 
nearest to the Sun along the line of sight, marked with P, because the source function of the 
scintillation signal drops down with the square of solar distance. In the case of helioglow 
observations, the maximum of the signal comes from the so-called Maximum Emissivity 
Region, which is located within 1.5-5 AU from the Sun, so it is important to carefully 
select the solar elongation of the line of sight. Note that the angular area of the polar region 
is quite small, so with an observations program that maps the entire sky only a small region 
in the maps indeed includes the signal from polar regions. 
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on the heliolatitude structure of solar wind are observations of IPS, interpreted by tomog- 
raphy modeling, that generally agree quite well with the Ulysses in situ measurements out 
of ecliptic and with the OMNI measurements in the ecliptic. Up to now, no continuous 
measurements of solar wind density as a function of ecliptic latitude have been available. 

Tests revealed that an appropriate balance between the latitudinal resolution of the cov- 
erage and fidelity of the results is obtained at a subdivision of the heliolatitudes' range into 
10-degree bins. Concerning the time resolution, the most welcome would be Carrington 
rotation averages, identical with the resolution of the photoionization rate and Lyman- 
alpha flux. Regrettably, such a high resolution seems to be hard to achieve because (1) the 
time coverage in the data from IPS has gaps that typically occur during ~ 4 months at the 
beginning of each year, which would induce an artificial 1-year periodicity in the data, and 
(2) the fast latitude scans by Ulysses were about 12 months long and hence differentiating 
between time and latitude effects in its measurements is challenging. Thus a reliable lati- 
tude structure of solar wind can only be obtained at a time scale of 1 year and this is the 
time resolution of the model we are going to present. 

Concerning the global mapping of solar wind parameters from ecliptic it has to be 
pointed out that the accuracy of measurements of solar wind parameters decreases with 
latitude because of geometry reasons. The polar values are the most uncertain (and possi- 
bly biased) because the signal in the polar lines of sight is only partly formed in the actual 
polar region of space, which can be understood from a sketch presented in Fig. 12.181 

In the following part of this chapter, we will present a procedure to retrieve the so- 
lar wind speed evolution in time and heliolatitude, and subsequently two procedures of 
retrieval of solar wind density. One is based on a correlation between the solar wind 
speed and density that was established from the three Ulysses fast latitude scans and must 
be regarded as an interim solution, to be used until the other one, based on the SWAN 
Lyman-alpha helioglow observations, will be available. 



2.5.5 Latitude profiles of solar wind velocity from Interplanetary Scin- 
tillation observations 

Interplanetary Scintillations (IPS) observations carried out by the Solar- Terrestrial En- 
vironment Laboratory (STEL) of Nagoya University (Japan) enable us to derive the lati- 
tude structure of solar wind speed and its variations in time. We use in the analysis data 
from 1990 to 2009, obtained from 3 antennas (Toyokawa, Fuji and Sugadaira), and from 
another antenna (Kiso) since 1994. The 4-antenna system w as operated until 2005, when 
the Toyokawa antenna was closed uokumaru et all 1201011 : since then the system again 
was operated in a 3-antenna setup. The IPS data from STEL are typically collected on 
a daily basis during 1 1 Carrington rotations per year: there is a break in winter because 
the antennas get covered with snow. The IPS observations are line of sight integration of 
solar wind speed weighted by density turbulances and Fresnel filter. Each day 30-40 lines 
of sight for selected scintillating radio sources are observed. The line of sight integration 
effect is deconvolved using the Computer Assisted Tomograph y (CAT) method developed 
by the STEL group iKoiima et al .1 120071: \Koiima et ali \\9m . The LOS's are projected 



on the source surface at 2.5 solar radii (R ), which is used as a reference surface in time 
sequence tomography. The latitude coverage of the sky by IPS observations is not uniform 
and is strongly correlated with the Sun position on sky, which changes during the year, 
and with the target distribution on sky, with relatively few of them near the solar poles. 
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Additionaly, the observations of the south pole are of lower quality than these of the north 
pole because of the low elevation of the Sun during winter in Japan. The original latitude 
coverage was improved owing to the new antenna added to the system in 1994 and by 
optimization of the choice of the targets. 

For our analysis we take solar wind speed from 1990 to 2009, mapped at the source 
surface on a grid of 1 1 x 360 x 180 records per year, which correspond to a series of Car- 
rington rotations. The data are organized in heliolatitude from 89.5 deg North to 89.5 deg 
South. 

A comparison of the tomography-derived solar wind speed with the in situ measure- 
ments by Ulysses we performed showed that the accuracy of the tomographic results de- 
pends on the number of IPS observations available for a given rotation. The intervals with 
a small number of data points clearly tend to underestimate the speed. Consequently, we 
removed from the data the Carrington rotations with the total number of points less than 
30 000. Small numbers of available observations typically happen at the beginning and at 
the end of the year and at the edges of data gaps. The selection of data by the total num- 
ber of points per rotation constrainted the data mainly to the summer and autumn months, 
when all latitudes are well covered. 

The selected subset of data was split into years, and within each yearly subset into 
19 heliolatitude bins, equally spaced from -90 deg to 90 deg. The speeds averaged over 
bins and over year for the latitudinal bins yield the yearly latitudinal profiles of solar wind 
speed, shown in Fig J2. 191 They cover a half of solar cycle 22 and full solar cycle 23. In the 
analysis we use the two-step calculation, first Carrington rotation averaged values per bin 
and next the yearly averaged calculated from the monthly ones. It is worth to note that the 
solar wind profiles for specific Carrington rotations have a very similar shape to the yearly 
profiles, which suggests that the latitude structure is stable during a year and changes only 
on a time scale comparable with solar activity variations. 
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Figure 2.19: Heliolatitude profiles of yearly averaged solar wind speed obtained from 
the Computer Assisted Tomography analysis of interplanetary scintillation observations 
0/1 1201 oh . 
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The results confirm that solar wind speed is bimodal during solar minimum, slow at 
latitudes close to solar equator (and thus the ecliptic plane) and fast at the poles. The lat- 
itude structure evolves with the solar acitivity cycle and becomes flatter when the activity 
is increasing. The structure is approximately homogeneous in latitude only during a short 
time interval during the peak of solar maximum, when at all latitudes solar wind is slow 
(see the panel for 2000 in Fig. 12.191 ). Shortly after the activity maximum the bimodal 
structure reappears and the fast wind at the poles is observed again, but switchovers from 
the slow to fast wind at the poles may still occur during the high activity period: compare 
the panels for 2001 and 2003 in the aforementioned figure and see the solar activity level 
depicted with the blue line in Fig. 12.161 During the descending and ascending phases of 
solar activity there is a wide band of slow solar wind that straddles the equator and ex- 
tends to midlatitudes; the fast wind is restricted to polar caps and upper midlatitudes. At 
solar minimum, the structure is sharp and stable during a few years straddling the turn of 
solar cycles, with high speed at the poles and at midlatitudes and a rapid decrease at the 
equator. Thus, apart from short time intervals at the maximum of solar activity, the solar 
wind structure close the poles is almost fiat, with a steady fast speed value typical for wide 
polar coronal holes, which is in perfect agreement w ith the measurements from Ulysses 
bPhillips etalih 995bl \McComas et all l2000bl l2006Tl . 

To verify the results obtained from the IPS analysis, we compared them with the data 
from the three Ulysses fast latitude scans and with the OMNI measurements in the ecliptic. 

The Ulysses velocity profiles used for this comparison were constructed from subsets 
of hourly averages available from NSSDC, split into identical heliolatitude bins as those 
used for the IPS data analysis and averaged. They are shown in Fig. 12.201 as red lines. 
Since the acquisition of the Ulysses profiles took one year each and the fast scans straddle 
the turn of the years, we show the IPS results for the years straddling the fast latitude scans; 
they are presented in blue and gray in Fig. 12.201 

The "grass-line" feature in the Ulysses profile for the second fast scan and in the equa- 
torial part of the first and third fast scans is due to the time Ulysses was taking to scan the 
10-degree bins. The fast scans were performed at the perihelion half of the Ulysses ellipti- 
cal orbit, with the perihelion close to the solar equator plane. Hence, the angular speed of 
its motion was highest close to equator and traversing the 10-degree bin took it less than 
one solar rotation period. Thus the "grass" is an effect of incomplete Carrington longitude 
coverage of the bimodal solar wind, with slow wind interleaved with fast wind streams. 
Near the poles the angular speed was slower and it took more than 1 Carrington rotation 
to scan the 10-deg bin and during solar maximum, when the gusty slow wind engulfed the 
whole space, the "grass effect" expanded into the full latitude span. By contrast, during 
the low-activity scans the solar wind speed at high latitude was stable, which resulted in 
the lack of the small-scale latitude variations in the Carrington rotation averages at high 
latitudes. 



2.5. Latitudinal structure of the solar wind and its evolution during the solar cycle 43 



800 
700 
600 
500 
400 
300 



800 



SW speed during 1st fast scan / 1994.7 - 1995.58 




-80 -60 -40 -20 20 40 60 80 
heliolatitude [degj 

SW speed during 2nd fast scan / 2000.9 - 2001.78 




-80 -60 -40 -20 20 40 60 
heliolatitude [degj 
SW speed during 3rd fast scan / 2007.1 - 2008.03 



800 




-60 - 40 - 20 20 40 60 
heliolatitude [degj 

Figure 2.20: Solar wind speed profiles from Ulysses measurements and IPS observations 
for the 3 fast Ulysses scans. Red: Ulysses, blue: IPS during the year of beginning of a 
Ulysses fast scan, gray: IPS during the year of the end of Ulysses fast scan. Top panel: 
the first fast scan during solar minimum, middel panel: the second fast scan during solar 
maximum, bottom panel: the third fast scan during minimum. 
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Figure 2.21: Yearly averages of solar wind speed from the IPS analysis (blue) and in-situ 
measurements collected in the OMNI-2 database (green), compared with the solar wind 
speed measured by Ulysses during its passage through the ecliptic plane during the three 
fast latitude scans (red dots). 



The profiles of solar wind speed from IPS and Ulysses are very similar, but some sys- 
tematic differences exist. On one hand, it seems that a ~ 50 km s _1 is a typical difference 
between Ulysses and IPS values in the polar regions, with the northern region usually in a 
better agreement than the southern. On the other hand, the agreement is sometimes almost 
perfect. 

The difference between the blue and gray lines in the top and middle panels of Fig. 
12.201 is a measure of true variation of the latitudinal profile of solar wind speed during 
one year. Ulysses was moving from south to north during the fast latitude scans, so the 
south limb of the profile from Ulysses ought to be closer to the south limb of the blue 
profile obtained from the IPS analysis, while the north limb of the Ulysses profile should 
agree better with the north limb of the gray IPS profile. Such a behavior is observed in 
the second panel, which corresponds to the solar maximum in 2001. In our opinion, it is 
a very interesting observation because it (1) shows how rapidly the latitude structure of 
solar wind varies during the maximum of solar activity, and (2) it confirms the credibility 
of both Ulysses and IPS results on solar wind structure. Originally, the interpretation 
of the speed profiles obtained from Ulysses was not clear, it was pondered whether the 
north-hemisphere increase in solar wind speed was a long-standing feature of the solar 
wind or was it just a time-variability of the wind at the north pole. Similarly, it was 
pondered whether IPS analysis is able to credibly reproduce the solar wind profiles given 
the fact that some of the profiles obtained approximately at the time of the fast scan seemed 
to disagree with the in situ data: our analysis suggests that the yearly-aveaged velocity 
profiles obtained from the tomography analysis of IPS observations agree with the in-situ 
observations from Ulysses even at times when the solar wind is restructuring rapidly at the 
peak of solar activity. 

The IPS data are in a very good agreement with the OMNI data collected in the ecliptic 
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plane (see Fig. l2.21t . Up to 2004 the agreement is almost perfect, afterwards small differ- 
ences appear. The agreement is better than with the in-cliptic Ulysses measurements from 
the fast latitude scans. This, in our opinion, is because the measurements of solar wind 
parameters at the ecliptic plane obtained during the fast latitude scans from Ulysses are 
challenging to directly compare with the OMNI and IPS measurements because Ulysses 
was passing ecliptic 10-degree latitude bin in a time equal to about half of the Carrington 
rotation and thus a reliable longitude averaging of the solar wind parameters could not be 
obtained. 

Thus the IPS solar wind speed profiles provide a solar wind latitude structure that can 
be adopted as an interim solution. They agree well both with the OMNI time series in the 
ecliptic and with the Ulysses measurements out of ecliptic for the time intervals when they 
can be compared directly. 
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Figure 2.22: Latitude profiles of yearly-averaged solar wind speed obtained from the in- 
terim procedure described in the text. The red lines show the speed profiles obtained from 
the Ulysses fast scans. 
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Figure 2.23: Correlation between the solar wind density and speed obtained from Ulysses 
fast latitude scans. Blue corresponds to the first scan (see Fig. 12.161 1. green to the third 
scan and red to the second scan, performed during solar maximum conditions. The dots 
correspond to the pairs of speeds and densities averaged over the 10-degree heliolatitude 
bins, the blue and green lines are the linear correlations specified in Eq. 12.401 The dotted 
red line is the density - speed relation proposed for the transition interval close to the solar 
minimum of 2002, calculated as the mean of the correlation relations obtained from the 
first and third latitude scans. 



For further analysis \Sokot et al.\ 1201211 smoothed the yearly speed profiles assuming 
that the variation of solar wind speed at the high latitude close to the poles is linear and the 
variation outside the polar caps can be approximated by a series of smoothly-transitioning 
parabolae. The resulting smoothed yearly heliolatitude profiles of solar wind speeds are 
shown in Fig. 12.221 These smoothed profiles will be used in the further part of this chapter 
to obtain the density and flux profiles of solar wind. 



2.5.6 Latitude structure of solar wind density and flux 

The time- and latitude-dependence of solar wind flux and density can be obtained using 
two different methods. The first one, which we regard as interim, relies on an approximate 
correlation between the solar wind speed and density inferred from the data from fast 
latitude scans by Ulysses. We discuss it first and show the results. The other one is 
expected to be the ultimate one and will be based on future results of analysis of the 
SWAN measurements of the Lyman-alpha helioglow, once it is completed. Here we will 
show just the method itself. 

Solar wind density and flux from density-speed correlation 

Solar wind speed and density seem to be correlated in heliolatitude, at least during the 
solar minimum conditions and for the observations collected during the fast latitude scans. 
The correlations seem to be slightly different between the first and third latitude scans, as 
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Solar wind density for 2. fast scan 
measured vs calculated 




-60 -40 -20 20 40 60 80 
heliolatitude [degj 



Figure 2.24: Comparison of solar wind density averaged over 10-degree heliolatitude bins 
actually measured by Ulysses (red) and calculated from the correlation formula in Eq. 12.401 
(gray). 



illustrated in Fig. 12.231 

"Ulysses O) = O-scan + b scan V (2.40) 

where a scan ,b scan are fitted separately for the speed and densit y values averaged over 10- 
degree bins using the ordinary least squares bisector method jlsobe et al. I ll 990ft (which 
allows for uncertainty in both ordinate and abscissa), separately for the first and third 
latitude scans. For the first scan (blue line and points in Fig. 12.231 ) we obtain a/,- rM = 
12.69, b first — —0.01332, for the third scan (green line and dots in Fig.[2 i 23]> the correlation 
formula parameters are a,^ = 10.01, b t urd — -0.01107. Thus, the slopes are almost 
identical and the main difference between the two formula is in the intercept, which reflects 
the overall secular decrease in solar wind density. 

The relation between density and speed for the second scan, which occurred during 
solar maximum, does not seem to be linear, but in this case the spatial and temporal effects 
seem to be convolved (as discussed earlier in this section). Therefore, we propose to use 
a relation that will be an arithmetic mean of the relations for the first and third scans: 
a-sec ond - [a first + athird) /2 and b second - (b f^st + b lhi rd) /2. This relation is shown in 
Fig. [2723] as the red broken line and a comparison of the density values actually measured 
during the second latitude scan and calculated from the correlation formula is shown in 

Fig. Em 

The interval of applicability of the latter formula is from 1998 until 2002. The formula 
from the first scan is applicable to the interval before 1998 and the formula from the third 
scan for the interval after 2002. 

We calculate the interim yearly profiles of solar wind density as a function of helio- 
latitude by applying Eq. ( 12.40b to the speed profiles obtained in the preceding subsection. 
Since this is an interim and very approximate solution, which by its nature is not very 
accurate, we do not attempt to further adjust it at the equator to the corresponding OMNI 
densities. The results are shown in Fig. 12.251 The accuracy of the results in the polar 
regions is on a level of 20% to 40%. 
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Figure 2.25: Latitude profiles of yearly averaged solar wind density obtained from the 
interim procedure. The red lines show the Ulysses fast scan profiles. 
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With the density and speed profiles on hand, one can easily calculate the flux: 



(2.41) 



dynamic pressure: 




(2.42) 



and charge exchange rate (in the stationary H atom approximation; Eq. 12.141 1. They are 
collectively shown in Fig. 12.261 

To obtain the solar wind parameters (speed, density, flux, dynamic pressure, charge 
exchange rate and total ionization rate) at a monthly resolution, we replace the equatorial 
bin directly with the Carrington rotation averaged series from OMNI, linearly interpolated 
to halves of Carrington rotations. The +10 deg bins are replaced with values linearly in- 
terpolated between the ±20 deg bins and the equatorial bin. The pole values are calculated 
from the the parabolic interpolation between the ±70 and ±80 deg bins, because due to the 
problems discussed earlier in this chapter direct measurements over the poles are not avail- 
able. The remaining latitudinal bins are linearly interpolated in time between the yearly 
profiles. As a result of such a treatment, we utilize all the available information on the 
equatorial bin of the solar wind. Away from the equatorial bin, where such an information 
is not available, we have a smooth transition into the region of low time-resolution model. 

Summing up this portion, we have the structure of solar wind speed from the smoothed 
IPS profiles (Fig. 12.221 1. The density structure is obtained (Fig. I2.251 l from the density- 
speed correlation from Ulysses (Eq. ( 12.401 ). Fig. 12.231 ). From these, we calculate the solar 
wind flux, dynamic pressure, and charge exchange rate between solar wind protons and 
neutral H atoms (Fig. 12.261 ). 
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Figure 2.26: Contour maps of solar wind flux (in cm 2 s ), charge exchange (in s ) and 
dynamic pressure (in nanopascals) shown as a function of time and heliolatitude. 
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Outlook: latitudinal structure of solar wind flux and density from IPS and Lyman- 
alpha helioglow observations (a sketch of the method) 

We assume that the total ionization rate from SWAN is the sum of the charge exchange 
and photoionization rates of H, with electron impact ionization neglected (because [3 S \ «: 
Pcx for the MER distance and beyond). Thus the heliolatitude- and time-dependent charge 
exchange rate can be calculated as 

Pcx (<t>j, ti) = Act (<t>j, ti) - Aiph (0, td (2.43) 

where yS to t (<pj, ti) is obtained from inversion of a SWAN map and the photoionization rate 
Aiph (0, f,), here assumed to be spherically symmetric, is obtained from one of the formulae 
specified in Eqs. l2.20ll2.21l or l2.22l 

Since the total ionization rate data from SWAN for the equator do not always agree 
with the rate derived from direct in situ measurements, we calculate the contrasts k^swan Wj> 
of the SWAN-derived ionization rates as a function of heliolatitude using: 

/ , /Wan yPjftn 
k/j,swan [4>j, ti) = -r— , (2.44) 

V PSWAN (U, ti) 

where jSswan (<Pj, ti) is an ionization rate at <^j-th heliolatitude and /Sswan (0, f,) is an ioniza- 
tion rate at the equator. In this way we can build yearly contrasts for the 19 heliolatitudes 
bins and for the years from 1996 to 2010. Now, multiplying the contrasts kb,swan (<Pj, ti) 
with the monthly averages of the total equatorial ionization rates (see Fig. 12.121 ). we obtain 
latitudinal profiles of the total ionization rate that agree with the baseline values for the 
equator. 

Now we can calculate the absolute charge exchange ionization rate for all heliolatitudes 
by subtracting the photoionization rate from the total rate Eq. (12.43b . In view of the formula 
for charge exchange rate given in Eq. ( 12.141 ), with the solar wind velocity profile obtained 
from IPS observations v (cpj, ti) we can now calculate the total solar wind flux as a function 
of heliolatitude in the following way: 

F (<f>j, tt) = yScx tt) /crcx (v (<pj, ti)) (2.45) 

With the flux and solar wind speed profile on hand, it is straightforward to calculate the 
profile of density: 

n(4>j,t,)=F(tj,ti)/v(<pj,ti) (2.46) 

In that way we can obtain a model of evolution of solar wind speed and density as a func- 
tion of time and heliolatitude that can be used to calculate the evolution of ionization rate 
for H atoms traveling at arbitrary speeds and to be input to global models of the helio- 
sphere to calculate the flux and dynamic pressure of solar wind. This calculation will be 
performed once the intercalibration of the heliospheric EUV measurements, which is the 
goal of this working team, is completed and the final inversion of the SWAN photometric 
observations for the ionization rate profiles is performed. 
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2.6 Summary 

Intercalibration of heliospheric UV and EUV measurements requires a common basis 
of heliospheric ionization processes, which affect the distribuition of neutral interstellar 
gas in the heliosphere and thus influence both the spectrum and intensity distribution of 
the heliospheric backscatter glow. In this chapter we present a review of the solar factors 
affecting the distribution of neutral interstellar gas in the heliosphere, namely we discuss 
the radiation pressure, solar EUV ionizing radiation, and solar wind. We review the history 
of measurements of these factors and develop a model of time and heliolatitude evolution 
of solar wind speed and density based on data available from in situ measurements of the 
solar wind parameters, remote sensing interplanetary scintillation observations of the solar 
wind speed structure, and from remote sensing observations of the heliospheric Lyman- 
alpha backscatter glow. The results of this model are used as input in the global helio- 
spheric models, discussed by Izmodenov et al. in another chapter of this book. The results 
of the global heliospheric modeling can in turn be used to fine tune the absolute calibration 
of heliospheric EUV measurements and, in a new iteration, to further refine the model of 
evolution of solar wind. The refined model will cover at least the two recent solar cy- 
cles and thus will provide a common homogeneous basis for interpretation of the present 
and past heliospheric experiments. It will also be used to interpret th e observations of 



Energe tic Neutral Atoms by the Interstellar Boundary Explorer (IBEX) \Mc Comas et al 
l2009bl latl. 
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